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ABSTRACT 



The general problem chosen for consideration concerns the 
establishment of a payload in orbit about the planet Mars Empha- 
sis in this study has been placed upon a comparison of the various 
paths available for a free-fall, space stabilized rocket, the re- 
quirements .necessary to achieve such paths, and the selection 
of the optimum path. The optimum path has been compared with 
other trajectories and conclusions have been drawn as to the 
requirements for deviations from the optimum path. The study 
is primarily concerned with the midcourse portion of the trip, 
which has been defined by the authors as that part of the overall 
trajectory extending from a desired circular orbit about Earth 
to a desired circular orbit about Mars. A firm basis for the 
selection of the orbit has been established. Treatment has also 
been given in this report to the co’ninuous determination of 
present position .in space and ’be initiation of corrective action 



when deviation from a desired path exists. The study also 
includes generalized instrumentation and guidance requirements. 
Extensive quantitative information has been compiled which will 
serve as a basis for the determination of instrumentation re- 
quirements as well as creata a proper background for an 
appreciation of flight within the solar system. 
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OBJECT 



The object of this paper is to investigate the mid-course 
guidance requirements to place a payload in orbit about Mars. 



CtL\PTER I 



INTRODUCTION 



The authors have undertaken this study in order to bring 
together, in proper context, some considerations which may be 
anticipated in interplanetary flight . A journey to the planet Mars 
has been chosen as the basis for the presentation. However, 
many of the ideas and methods set forth here can be considered 
applicable to flight to other bodies within the solar system. It 
is felt that no loss of generality will result if attention is focused 
upon this particular choice. 

This paper has evolved as a study of several problems as- 
sociated with interplanetary flight. The objective to be pursued 
is concerned with placing a payload in orbit about the planet Mars . 
Departure is assumed to originate from orbit about Earth. The 
attainment or possibility of flight is not questioned. Rather, 
this is a broad study of some preferred methods, minimum re- 
quirements, and ideas concerned with fulfillment of the general 
objective . 

Recent developments have justified studies of interplanetary 
flight. It is realized, however that the topic is hopefully opti- 
mistic in many respects. For this reason the study is neces- 
sarily broad and is not designed to give an intricate, detailed 
account of methodology 

A flight to the moon will logically precede planetary flight, 
but its close proximity causes this problem to become quite 
specialized as contrasted to the Mars flight. Piloting a boat in 
an inland waterway as compared to navigating a ship on the open 
sea is an appropriate analogy. The two problems are independent 



in a large sense . 

« 

This study has been preceded by an extensive search through 
the literature for information regarding this comparatively new 
field. Treatments have generally varied from the popularized 
texts to the highly specialized mathematical devdopments. 

There appears to be a paucity of practical approaches to the 
problem which can be interpreted by the systems designer or 
engineer. 

The authors generally envisage an unmanned vehicle with 
small payload, unpowered except for terminal and corrective 
thrust requirements . It is believed instrumentation can be 
handled with components whose performances do not differ too 
radically from slight extensions of present day development. 

A practical simplified view has been taken as a rule due 
to the lack of an available high speed digital computer. Such 
a device would of necessity be required if accurate trajectory 
studies were in order. However, it is felt that little departure 
from reality exists under the assumptions imposed. In fact, 
the simplifications add to the comprehension of the material. 

Solutions are sought in answer to the following questions: 

(1) What are the major influences of the solar system 
upon the rocket in flight ? 

(2) What simplifying assumptions can be made? 

(3) What routes may be taken? 

(4) lYhat criteria are available for the selection of 
one trajectory or route in preference to another? 

(5) How may position in space be continuously 
determined for navigational purposes? 

(6) How may the rocket be assured of arriving in the 
target area if errors do occur? 
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CHAPTER 2 



THE SOLAR SYSTEM 



A summary of currently accepted data concerning the 
universe, which is pertinant to this study, is needed for formu- 
lation of assumptions and as basis for computation. This is 

presented below as discussion and in tabular form and figures, 

2 4 Q 44 ‘=>4 'iQ* 

and has been derived from references. ' ' ' ’ ' ' 

2. 1 The Sun and its Planets 

The solar system consists of the sun, nine conspicuous 
planets, and numerous smaller planets. The conspicuous planets 
are placed in two groups. The major planets , Jupiter, Saturn, 
Uranus, and Neptune, are large, have low density, and are 
surrounded by high atmospheres. The terrestrial or minor 
planets . Mercury, Venus, Earth, Mars, and Pluto are similar 
to the earth in many respects. The so called navigational planets 
are Venus, Mars, Jupiter, and Saturn, deriving this classi- 

9 

fication from their nautical use for position determination. 

For the period of 1956-1970, the angular motion of the 
major planets in the solar system will be as shown in Fi^. 2.1 . 
The minor planets. Mars, Earth, and Venus, are shown in their 
positions as of 1 January 1970, and in their respective orbits, 
drawn as circles. 

A generalization which serves to illustrate the area occupied 
by the planets and sun is to picture a piece of 1/4 inch plywood cut 
in a circle of 5 feet diameter. This would contain all the planets 
in their orbits around the sun (except Pluto) if drawn to a scale 
such that Neptune's major planetary axis was 5 feet. 

* Superscript numbers refer to bibliography. 
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Fig. 2.1 Angular motion of major planets from 1965 to 1970 
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2. 2 Solar Radiation and Temperature in Space 



The solar constant of sun’s radiation is 1. 938 calories per 

54 

minute per square centimeter, ±.05 This value is equiva- 
0 

lent to 1.35 X 10 ergs per second per centimeter square or 
0 

87.1 X 10 ergs per second per inch square at a distance of 1 AU 
from the sun. This incident energy varies as the reciprocal of 
the range squared, as shown in Fig. 2.2 

The maximum temperature of the surface of a black body 
planet with no atmosphere varies as the reciprocal of the square 
root of the distance of the body from the sun, as shown in fig. 2. 3. 
At a distance of 1 AU from the sun this temperature is 392°K. 

2. 3 Approximate Orbital Paths of Planets 

The path of a planet or body in orbit about the sun can be 
completely described by seven elements, in so far as the body 
depends only upon the sun's attractive force. These are: 

1 . Semi ma.jor axis (a) 

2. Eccentricity (e) 

3. Inclination (i) 

4. Longitude of ascending node 

5 . Longitude of perihelion ( n) 

6 . Period (T ) 

7 . Epoch 

These are shown in Fig . 2.4 and 2.5. 

Kepler's Three Laws of Planetary Motion are also appli- 
cable and are: 

1. The planets move in ellipses, with sun at one of 
the focuses 

2 , The radius vector sweeps out equal areas in equal 
time 
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Fig. 2.2 Incident solar energy as a function of distance from Sun 
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Black body temperature as a function of distance from Sun 




Fig. 2.4 Basic ellipse 




Fig. 2.5 Basic ellipse oriented to reference plane 
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3. The square of the time of revolution is proportional 
to the cube of the mean distance. 

2.4 Perturbations of Planetary Orbits 

The actual motion of the planets may best be described as 
follows; 

1. The orbits of the planets are not fixed but gradually 
change in eccentricity, inclination, (etc. ). These 
changes are classified as secular p erturbations ; 
however, these perturbations are so slow that during 
any one revolution the orbit is altered very little and 
during the time of any given single interplanetary 
flight, can be neglected . 

2 . A planet in its motion does not exactly follow 
Kepler's Laws but is seen to vary in three dimensions 
relative to its ideal Keplerian Ellipse. Such changes 
are very small and are called periodic perturbation. 
These values expressed in heliocentric longitude (Xp) 
are: 



Mercury 


15" 


Venus 


30" 


Earth 


1’ 


Mars 


2' 


Jupiter 


30' 


Saturn 


70' 


Uranus 


60' 


Neptune 


35' 



In general the larger the terms the longer the period of 
periodic perturbation. For example Jupiter and Saturn exhibit 
periods of periodic perturbations of about 913 years . 
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2. 5 The Invariable Plane 



The most stable property of the solar system is the so 

called invariable plane. This is defined by the condition that the 

— — — — ^ 

total angular momentum of the solar system is zero about any 
axis in this plane. The invariable plane can be considered as 
the fundamental reference plane of the solar system and is 
inclined 1*^ 38' 7" to the ecliptic reference plane of 1900 with 
ascending node in longitude 106° 42' measured fromT* . This 
inclination is decreasing at the rate of .144" per year and will 
reach 47' in about 20, 000 years at which time it will begin to 
increase again. 

2 . 6 Cislunar Motion 

Cislunar (earth -moon) motion in the plane of the ecliptic is 
that of two unequal revolving masses in a central force field, 
whose center of mass follows a slightly distorted Keplerian 
ellipse of eccentricity (e)=. 01674. This cislunar center of mass is 
approximately 2880 miles from the earth center in the direction 
of the moon. See Fig. 2. 6 and 2.7. 

The moonte path around the earth has an eccentricity .055 
and a mean inclination to the ecliptic of 5° 9'. The average 
distance from earth to moon is 238^. 857 miles (60.207 earth radii) 
with maximum distance of 252, 710 miles and minimum distance of 
221, 463 miles . 

If we represent the cislunar orbit by a circle of radius 100 
feet, the moon will deviate from this circle approximately 3 inches 
on either side in one transit around the earth. During this time 
the geocenter will deviate approximately .04 inches either side 
of the orbit. 

2. 7 Solar Motion 

Campbell and Moore in 1920 found the solar system to be 
moving at about 12.2 miles per second wilh respect to the nearer 
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2. 5 The Invariable Plane 
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moving at about 12.2 miles per second with respect to the nearer 




(T) 

FIRST POINT OF ARIES 

/PRECESSES 50.26l97yr\ 

\T= 26,000yrs / 

Fig. 2, 6 Motion of theEarth in the ecliptic plane 




Fig. 2. 7 Cislunar Motion 
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few hundred stars toward the constellation Hercules , Jay, 

Sanford, Stromberg, and Hubble, among others, estimate that 
the solar system and the neighboring stars are moving at about 
170 miles per second toward the constellation Cygnus , 

Conclusions reached by Oort place the solar system about 
32, 000 light years from the center of the galaxy which is in 
Sagittarius and that this position is about two thirds of the dis- 
tance from center to circumference. At 170 miles per second 
the period of galactic revolution of our sun is approximately 
2, 200, 000 years, indicating that during the time covered by any 
given flight within the solar system, the sun may be considered 
as a rather stable reference point. 

2 . 8 Fixed Stars 

The nearest star. Alpha (,’entauri, is at a distance from our 

13 

sun of 4.3 light years or 2. 53 x 10 miles. It has also been de- 
termined that no star exhibits a parallax of more than . 76 seconds 
and only 200 stars have proper motion greater than one tenth of 
"Barnard's Runaway Star" which traverses 10.25 seconds of arc 
per year. For the purposes of this paper a "Fixed" star, as the 
term is used for establishment of celestial inertial reference, 
will be considered to be a star exhibiting a parallax and a proper 
motion of less than one quarter second of arc.Ofthe stars checked 
to date only about . 01% do not fall within this Classification of 
"Fixed" stars, 

2. 9 Gravitational Effects 

Newton's general Law of Gravitation states that, "Any two 
particles of matter attract each other with a force proportional 
to the product of their masses and inversly proportional to the 
square of the distance between them"^"^ or f = K where K is 

the so called Gaussian gravitational constant equal to 



6.67 X 10 



. From this the acceleration of gravity 



-8 



3 



cm 




KM 



is then 



g = 



2 • 



R 



Values of the acceleration of gravity for the various solar 
bodies are tabulated in Table 2. 2 with the mean surface 
gravity for earth equal to 1, Fig. 2. 8 presents "g" versus the 
radial distance from the earth, mon and Mars and the variation 
in the effect of the sun's gravitational attraction over the distance 
between the orbits of Earth and Mars . 

18 

Herrick has pointed out that Gauss was able to determine 

2 

mathematically the value of K , for heliocentric orbits with dis- 
tance in AU, accurate to nine significant figures but neither the 
laboratory value of gravitational acceleration nor laboratory units 
of measurements at present will permit any greater accuracy than 
3 or 4 significant figures. 

An interesting discussion of Gravitational and tidal effects 

which has direct bearing on any travel between Earth and Mars 

56 

has been presented by Schaub . Schaub states that the inner 

54 

most moon of Mars, Phobos, is within RockeLs limit . This 
means that since Phobos is assumed to have the same density as 
Mars and is traveling in an orbit of less than 2. 44 times Mars 
radius (see Table 2.2 ) its surface is now free of all loose 

material. Eventually, due to the bombardment by rays from 
space, more and more of this moon will be torn away forming a 
ring of dust particles around Mars . This implies that any 
terminal orbit about Mars should stay well clear of Phobos less it 
hasten this disintegration possibly destroying the vehicle in the 
process . 
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Fig. 2. 8 Gravitational field effects 
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CHAPTER 3 



THE PROBLEM SIMPLIFIED 



3.1 Introduction 

The Law of Universal Gravitation as enunciated by Newton 
states "Every particle of matter in the universe attracts every 
other particle with a force which acts in a line joining them, and 
whose intensity varies as the product of the masses and inversely 
as the squares of their distances apart. " 

44 

Moulton presents an excellent discussion of this law and 
states, "it will be observed that the law of gravitation involves 
considerably more than can be derived from Kepler's Laws of 
Motion; and it was by a master stroke of genius that Newton 
grasped it in its immense generality and stated it so exactly that 
it has not been necessary to change a single syllable in over two 
hundred years. When contemplated in its entirety it is one of the 
grandest conceptions in the physical sciences." 

It is evident that the various bodies comprising the Solar 
System will exert certain forces upon a rocket in flight among 
them. Because of the extremely large masses and distances 
involved we will omit consideration of the minor bodies, such as 
asteroids and meteor swarms, and show later by inference that 
this omission is justifiable. Considering then only the Sun, Moon, 
rocket, and nine planets, we are faced with a twelve -body problem. 
It would appear expedient to simplify this problem if possible 
since twelve bodies require thirty-six second order differential 
equations (thirty -three, if referred to the Sun) to adequately de- 
scribe their motions. 
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3. 2 Differential Equation of Motion of the Rocket 



In order to reduce the complexity of the problem let us first 

examine the terms of the equation of motion of the rocket and 

consider their magnitudes. With reference tOvFig. 3.1 define 

71 

I as an inertial reference frame in the convention of Wrigley 
and 0 as a frame within the Solar System non-rotating with respect 
to I and with axes parallel to those of I. Note that at this point 0 
is restricted only to lie somewhere in the vicinity of the Solar 
System. A represents the position of the rocket^ a variable. 

We may write Newton's Second Law in vector notation as 
follows; 



[R 



lAM 



1 = ^ 

Jt 



^f ^nf 



(3-1) 



where [ Rj^] j is the second derivative of the vector with 
respect to the frame I and the other terms are previously defined. 



Then; 



f 



f 



^0A ^CA + ^$A ®9A ^cTa 

+ + G^^ + G^^ + Gg^ (3-2) 



where G. . 

3 A 

If there 



is defined in Derivation Summaryl . 
is a particle at point 0 of mass mj then; 



• « 

m [ Rj ] |- = m (GL + GL + G^ + G^ 
o’- lo^ I o H0O ®o ao Qo 

and the masses may be cancelled. 

By geometry 



) 

(3-3) 



f ^lA^ I ~ t a] T 



oAJ I 



ToJi 



(3-4) 
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Fig. 3.2 Planetary configurations for Table 3.1 
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Hence by substitution 



[R 



lo^ I 






= ff + 



nf 



(3-5) 



or 



(G ^ + G^ + G^ 

' ©o ©o Co 






+ ^©A + ^dA 



)+[R Aj = ff 
^ o A' I nf 

. . . . ) (3-6) 



Hence in full 



f I " ^nf + ^^©A “ + ^^©A " + ^^SA '^do^ 

+(G2,a - G^jq) + (Gp^^ - Gy^o)+ (G^^- ^q) 

+(G^^^ - G^o^ + ^^0A ■ (3-7) 



We shall now consider applying this derivation (frame 0) to 
various members of the Solar System. 



3. 3 Inertial Points Within the Solar System 

Looking first at the motions of the Sun, or the Solar System 
as a unit, it appears from a statistical analysis of the motions 
of the "nearer" stars that the Solar System possesses a velocity 
of approximately 12.2 miles per second relative to these "nearer" 
stars. The "nearer" stars range from 4 to 300 light-years 
distant from the Sun, Further, the Sun is thought to occupy a 
position approximately 20, 000 light-years from the center of our 
Galaxy (The Milky Way) with a period of revolution about the 
galactic center on the order of 2, 200, 000 years. (See Section 2. 7 
for further details on the Solar System) . In addition our Galaxy 
is a member of the so-called Local Group consisting of 17 galax- 
ies, clusters, and clouds including the Andromeda Galaxy, The 



34 



Large and Small Magellanic Clouds, and several smadl clusters. 

It is thought that these 17 members rotate about a common center 
of mass located somewhere between our galaxy and the Andromeda 
Galaxy. Lastly, current theories on the "Red Shift" of the optical 
spectra of distant galaxies conclude that these galaxies are re- 
ceding from the Solar System at such tremendous relative ve- 
locities that certain effects attributed to the Theory of Relativity 
must surely be present. It is not the purpose of this paper to 
argue any of the current theories or to attempt any analysis in- 
cluding Relativity Effects, These motions to which the Solar 
System is subject are mentioned solely as a basis for the following 
statements: 

The Sun is probably undergoing motions involving intricate and 
complex accelerations; however, the vast macrocosmic scale of 
distances involved is such that during the time of flight of our 
rocket the effect of these ill-defined accelerations is neither ap- 
preciable nor measurable. This then permits the statement, 
"Newtonian space is that unaccelerated with respect to the 'fixed' 
stars", particularly since it is not known that any stars are 
fixed, but only that their motions are very small during recorded 
history. 

The foregoing definition of Newtonian space is preferred by 
the authors since it appears at least semantically more defi- 
nitive than the statement "Newtonian space is that space wherein 
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Newton's Laws are valid." We say then "The possible motions 
of the center of the Solar System, however complex, within our 
detective ability are not known to interfere measurably with the 
information obtained while ignoring these motions and effects 
thereof. " 

3. 4 Simplified Differential Equations of Rocket Motion 

On the foregoing basis we allow ourselves to define the center 
of (1) the Sun, (2) the Earth, and (3) Mars each in turn as the lo- 



cation of our inertial reference frame 0. Equation (3 t 7) may 
then be written successively; 



" ^nf ‘ ^ ‘ ‘ ox 

( 3 " 0 ) 

f ^©A^ I " \f ^©A ■’■^^OA" %©^ ^^(LA ~ ^(L©^ ' ' '(o Q\ 

\ 3- y ) 

[ ^dA^ I " ^nf ^CfA'^ ^^OA " ^OCf) + - G(^^) . . . 

( 3 - 10 ) 

2 

Table 3.1 gives values in ft/sec of the maximum values 
the above terms in parenthesis will hold, using mean values for 
the orbital radii of the bodies. Use of these mean radii is con- 
sidered justifiable since only rarely will the bodies be aligned 
such that these maximum values actually occur. Figure 3. 2 
clarifies the manner in which these values do occur. 

By inspection of Table 3.1; if we establish that all acceler- 
— 6 2 

ations less than 10 ft /sec are to be ignored, we may eliminate 

from future consideration most of the parenthetic terms in the 
foregoing equations. We then have; 

©zJ I " ^nf ^©A ■^^^©A ^^(LA “ ^€©^ 

+ (3-11) 

^ ^©A^ I ^ ^nf ^®A '''^^©A " ^©©^ ^ 

"^©^ (3-12) 

I " ^nf %A ^^©A" A" ^<rcf^ 

'*’^^©A~ (3-13) 

We have now reduced the problem to one of five bodies rather 

than twelve. We will investigate the possibility of reducing it even 
further. 
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R4_4| = 753, 000 " 

Fig. 3. 3 Illustration of Rocket -Moon Separation 
3. 5 The Lunar Gravity Field Effect 
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Referring to the Hohmann Ellipse it may be shown that in 

the first 15.2 days of rocket flight along this trajectory the radius 

vector of the rocket measured from the center of the Earth 

increases by approximately 600, 000 miles. During this period 

the rocket is in almost linear flight, hence the magnitude of the 

Lunar gravity field is easily estimated and is shown in Fig. 3.3. 

At the end of the 15. 2 days the value of the gravity field has de- 

-5 2 

creased to less than 1 x 10 ft/sec and is continually decreasing. 

As a result of these rudimentary estimations it is evident 
that the effect of the lunar gravity field need be considered only 
during the initial portion of the trajectory. The example selected 
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above is unique in that it has been shown to require the least possi- 
ble rocket velocity hence the rocket would be subjected to Lunar 
influence for a longer time than in any other trajectory. The fore- 
going estimates then are an upper bound to the values which actu- 
ally may occur . v 

3. 6 The Flight-Phase Concept 

Similar arguments to those of the preceding section may be 
employed to justify ignoring the effect of the Earth and Mars grav- 
ity fields over certain portions of the flight. Certainly if a com- 
puter is designed to employ the values of these fields there will 
be a certain sensitivity level below which no cognizance of these 
values will be taken. 

Thus we may write with some degree of justification the 
following equations; 



f I " ^nf 




^©A 


Phase II 


(3-14) 


f I ^ ^nf 


-h 


^©A 


Phase I 


(3-15) 


t I " ^nf 




^d'A 


Phase III 


(3-16) 



where the limits of the Phases are to be identified. 



Thus by consideration of the trajectory and the magnitudes 
of the forces involved we reduce our problem from one of 12 bodies 
to three separate two -body problems. 



TABLE 3.1 



This table displays values for the parenthetic terms of 
equation 3- 7 In the general form we have; 



|(G. G. )|= k. R .^C — 
I jA 30 I 3 03 ^ 
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R. 

JO 



where k^. is the surface gravity of the indicated body 
R . is the radius of the body in AU 

Rj^ is the distance separating the body and the rocket in AU 

R. is the distance separating the body and point 0 in AU 
J o 
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0 = ® 


0 = cf 


|g©a-Go I 
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2.78 X 10"^ 
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-4 
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1 05^-09 1 
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tl.Ol X 10"^^ 


tl.Ol X 10"^^ 



* by judicious choice of firing time 
t Assuming that Pluto and Earth have the same 
value for surface gravity, 
t Sample calculations are given in Table 3. 2. 



SP 



TABLE 3 2 



Sample Calculation for the Values Listed in Table 3 1 
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CHAPTER 4 



THE TRAJECTORY 
A. SOME POSSIBLE TRAJECTORIES 

4..1 Basic Criterion for Selection 

Selection of the type of trajectory to be followed during the 
mid-course phase of an interplanetary exploratory trip of any 
type is dependent upon the basic criterion of when, in the future 
we contemplate attempting the flight. K we desire to make the 
attempt within less than ten years it is reasonable to postulate 
that presently developeci power plants and guidance systems must 
be used with considerable refinement a necessity. If the attempt 
will not be made until some time after 1970, more exotic power 
plants and guidance can probably be postulated. However, the 
basic trajectories can be summarized and a selection made for 
further study in this paper. 

4. 2 Classifications of Trajectories 

All trajectories, based on thrust available, can be divided 
into three classifications: (1) Powered, (2) Unpowered and (3) 
Combination of Powered and Unpowered. Although this statement 
is simple, the actual engineering and instrumentation re- 
quirements are exceedingly complex and somewhat dissimilar. 

4. 3 Powered Traje ctories 

The term "Powered Trajectory" implies continuous thrust with 
either acceleration or deceleration used during the entire flight. 
Also implied is the ability to apply this thrust in any direction by 
multiple thrust units or vehicle reorientation with respect to the 
trajectory. 
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The upper limit of this type trajectory can be considered to 

be strictly interceptor, while the lower limit is that of a distorted 

Keplerian Ellipse. Interceptor trajectories, using proportion- 
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al navigation, collision course, line of sight, or constant lead 
similar to present day antiaircraft and air-to^ir guided missiles^ 
require nearly unlimited thrust. However, many simplifications 
in guidance methods can be realized due to reduced time of flight 
and reduction in accuracy requirements if high speed inter- 
ceptor methods are used. 

With small thrust power units such as the ion gun, solar 

4T G 2 

furnace, or Nuclear-chemical, etc, ^ the trajectory can be 

14 

a logarithmic spiral such as developed by Forbes or a series 
of constantly changing ellipses with the actual trajectory being 
the envelope of the series. In addition the vehicle could be 
forced to traverse a predetermined elliptical path. This latter 
method would be extremely wasteful of both fuel and time . 

Powered trajectories are difficult to compute and almost 
infinite in number. The practicality of detailed powered 
trajectory study with any more than a general overall look 
at the field appears to be limited at present until the perfection of 
suitable thrust devices warrant the effort. 

4. i Unpowered Trajectories 

Basically the unpowered trajectory is a ballistic flight - 

22 43 72 

incorporating the Fire Control Problem ' with that of 

21 ^-4 60 

determining the path of a comet in space. ' " ' This implies 

the application of a large impulsive thrust applied over a very 

short interval of time with respect to total time of flight. This 

would result in a hyperbolic path with respect to a planet during 

/ 1.7 6 g 

application of thrust. " " At the end of this time the hyperbolic 
trajectory would merge with an elliptical path with respect to the 
sun, and the vehicle would be in free flight, acted upon basically 
by gravitational attraction of the sun with slight perturbations 
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caused by the planets, (see Section 3.6). 

A trajectory of this sort leads to simplification or 
elimination of a guidance system for mid-=course. However, 
extreme accuracy in the orientation of the velocity vector at 
termination of impulse is mandatory. Very small errors in 
magnitude or direction, or unforseen gravitational gradients 
during transit can cause large errors at the trajectory terminus 
as shown in Section. 8.6. 

4.5 Comb ination Trajectories 

Trajectories containing thrust periods and ballistic no- 
thrust periods are complicated by virtue of attempting to combine 
two Separately difficult problems into one. This method does 
permit slightly lower accuracies of the initial velocity vector 
but requires accurate mid=course position determination to be 
of any use. 

The basic advantage of this method is that presently developed 

thrust units can be used for correcting or changing the trajectory 

during mid-course, thereby increasing the probability that the 

mission will be completed. Using contemplated low thrust cor- 
47 62 

rective units ' powered phases could be of long duration and 
be equally as effective as high thrust, short burning time chemical 
units . 

4, 6 Trajectory Selected for Fur ther Study 

The basic mid-course trajectory chosen for further study in 
this paper is that of a combination of unpowered and powered. 

The ideal or hoped for situation is one in which no corrective 
thrust is needed during the time of transit. This would mean that 
a relatively unperturbed ellipse could be flown within the limits of 
terminal accuracy requirements. However, if it becomes ap- 
parent during transit that corrective action is needed, a new 
ellipse is computed and proper thrust applied to get on the new 
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trajectory. The position determination accuracy requirements 
for such a system are rather high as are the thrust control 
requirements. However, it is felt that these requirements are 
not entirely beyond present day engineering capabilities. 
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B, THE UNPOWERED TRAJECTORY 



4. 7 latroduction 

To intelligently decide upon a trajectory to be traveled to Mars 
one must weigh many aspects of the problem. In previous por- 
tions of this paper we have investigated some properties which 
influence a vehicle in flight in the Solar System; and, in addition, 
investigated some vehicular behaviors and parameters which 
could be selected or adjusted to our benefit. 

As previously described, the "unpowered" trajectory is a 
trajectory provided with an initial impulsive velocity, which, when 
perfectly instituted, needs no corrective action. This type tra- 
jectory lends itself well to analysis since the vehicle is essentially 
a body acting in a central force field (subject to perturbations). 
Primarily due to considerations of available fuel and performance 
characteristics it was decided early in this work to proceed with 
an analysis of this type trajectory. It is emphasized here that 
this type trajectory may by no means be considered more at- 
tractive than any other because of any inherent property associ- 
ated with it. It simply appears more feasible due to the current 
state of the art in the development of rocket vehicles. 
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Hohmainn appears to be the first author (circa 1925) to 
seriously consider the solution of space flight problems. One of 
the most important results of his excellent study is the description 
of the minimum energy ellipse, which as the name implies re- 
quires the least rocket velocity of any configuration of two planets. 
This trajectory has become known as the "Hohmann" ellipse and 
has the characteristic that the arrival of the rocket at Mars 
occurs when Mars is diametrically opposite Earth position at 



time of firing in their respective orbits about the Sun. Note 
Fig. 4.1 . 
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Preston -Thomas ' agrees with Hohmann and shows some 
simple relations of the Hohmann ellipse with other ellipses. 

Other than this the literature shows a decided lack of information 
on other trajectories of the "unpowered" type. 

It is the purpose of this chapter to show some relations exist- 
ing among the "unpowered" trajectories. These are the tra- 
jectories which when corrected for errors are described as the 
"combination" trajectories of Section 4,5 



4. 8 Basic Relations 

The basic relations of a vehicle in motion under the influence 

44 60 

of a central force field are set forth by Moulton , Smart , 
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Kooy and Uytenbogaart , and others. Those general equations 
selected by the authors which appear to best fit the problem are: 



r= a(l-e^) 
1 -e cos 6 
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e sin 6 
1 -e cos 6 



A fourth equation 



which will be of interest is 



evl-e^ sin 6 
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A general trajectory situation is shown in plan view in Fig. 4. 2 
Figs. 4 . 2 and 4 ^ 3 illustrate the quantities appearing in the 
foregoing equations . 



9 Some Simplified Trajectory Parameters 

We logically begin our analysis by considering the effect of 
the most massive member of the Solar System, the Sun, Re- 
ferring to Fig . 4 2 ; for a given trajectory we may establish 
initial and final positional conditions (i.e, - we know we are on 
the Earth and hopefully we will arrive at Mars). Hence we write 




Hohmann Ellipse 



Fig. 4.1 The Hohmann Ellipse 



(4-1) as 



r 




(4-5) 



o 1 - e cos 6 



o 



and 



r 



a(l - e^) 



(4-6) 



A 1 - e cos 6 



A 
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rocket ellipse 
aphelion 



Fig. 4.2 General t rajectory configuration 




Fig. 4.3 Some vector-angle relations for the general 
trajectory 



48 



where = | Rq^ | at t = t^, and =| Rocjjat t = 
Similar definitions apply to the angles 0^ and 0^. 



Let r =1 AU 
o 



Let r^ = 1.5237 AU 



Dividing (4-5) by (4-6) and rearranging, we have 



e = 



0.5237 



1. 5237 cos 0 A = cos 0 
A o 



(4-7) 



Then from (4-1) 



a = 



(1 - e cos 0 ) r 
o o 



1 - e 



(4-8) 



which may be expressed more conveniently 
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1 2 
1 - e 


a 


1 - e cos 0 

o 



(4-9) 



where a is in AU . 

This defines two of the parameters (a and e) of any ellipse. 

We may novr establish the vector velocity of the rocket with 
respect to the Sun in terms of the trajectory parameters. 






at t = t^. Then (4-2) may be written 



V = 97, 500 ^/2 - 1/a 



ft/sec 



where a is in AU, 



(4-10) 



Let be the angle shown in Fig. 4.4 at t = t^. Then 



<t> = tan 

^o 



-1 



e sin 0 



o 



1 “ e cos 0 



o 



(4-11) 
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Fig. 4.5 Velocities relative to Mars 
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Additionally, from Kepler's Laws 




(4-12) 



and the period of the ellipse is known. 

The boxed equations summarize pertinent parameters in the 
study of rocket motion with respect to the Sun, and form a basis 
for the calculations of velocity vector requirements with respect 
to Earth and Mars . 

It might be argued here that greater accuracy should have 
been demanded of these basic equations (i.e. - Use 1.523688 AU 



Known accurately only to about one part in 10 hence no good 

•% 

purpose could be served by extending the required accuracy of 
these equations beyond the fourth decimal, at least for a pre- 
liminary analysis . 

Calculations were carried out for about 50 elliptical tra- 
jectories and the values for e andj j are plotted in Fig 4. 6 
and 4. 7 as aids in future studies. These plots presuppose that 
6 and 0 . are given, known, or selected. A logical basis for 
the selection of 6 ^ and 0^ is the subject of another section. 
Ellipticity < 1 was arbitrarily chosen since excessively high 
rocket velocities are required for parabolic or hyperbolic tra- 
jectories . 

4 10 Rocket Velocity Vectors with respect to Earth and Mars 

In the preceding section equations were developed to express 
parameters of the elliptical trajectory and initial velocity of the 
rocket with respect to the Sun. Obviously, since Earth and Mars 
are also bodies in motion with respect to the Sun we want the 
velocity of the rocket with respect to these bodies. The relations 
are simple, geometric in nature, and are developed herewith. 




The astronomical unit is 
. . ,«4 , 



At t 



o 
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Fig. 4.6 Ellipticity as a function of 6 and 6 
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Fig. 4.7 Rocket velocity magnitude with respect to the Sun as 

function of and 6 . 



= V0A + Vq^ 



(4-13) 



At t = t + tj, 
o f 






as defined in Definition Summary 1 
From equation (4-3) 



d> = tan 
^o 



-1 



e sin 0 



1 - e cos 0 



o 



And from Fig. 4.4 



't . -1 



V sin 6 

o 



V cos 6 - Voi 

o ^o tp 



Vr 



V cos 6 - Vr\oi 



0A ■ cos 11/0^ 

The properties of an ellipse identify 



ifV-A= Vqa 



Va = 


r 

o 


'V 


A 


’^A 


o 


at t 


= t 





<l> j^= tan 



-1 



o [ 
e sin 0 , 



1 - e cos 0 



AJ 



and from Fig. 4. 5 and equation (4-14). 






tan 



-1 






'^A 



(4-14) 



(4-15) 



(4-16) 



(4-17) 



(4-18) 



(4-19) 



(4-20) 
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VC^A- 



(4-21) 



cosx|;-C^^ 

If r = 1 . QOOO AU and r * = 1 . 5237 AU then 
o A 

Vq© = 97, 500 ft/sec 

and 

990 ft/sec 

Figures 4. 8 and 4.9 show ^ and as functions of 0^ 

and 0^. These values are important primarily as a step in 
solving the trajectory problem and are given here for possible 
future use by the reader. 



4. 11 Time of Flight 



The calculation of time of flight is necessary at this point in 
order to determine the angular separation of Earth and Mars at 
the instant of rocket firing. Furthermore, it is necessary to 
know tj in sidereal or calendar time since the inclination of Mar’s 
orbit with the ecliptic is a function of real time. 



If we substitute t^ and t^ as defined in Definition Summary 1 

each in turn in Equation (4-4) and recognize 

0=0, we haVe 
o d 



J, 2 . 

Wl-e SI 



sin 0 



1 - e cos 0 

o 



sin"^ 



e - cos 0 
£ 

- e cos 0 



o 



^2 




(4-22) 



and 



iJl- 



e sin 0 



1 - e cos 0 



A ^ . ~1 

+ sin 



A 



e - cos 0^ 

1 - e cos 0, 



loMrs 
+l-=t ■/ 



(4-23) 



Substituting tf = after subtracting (4-23) from(4- 22) we 

have 
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Fig. 4. 8 Rocket velocity magnitude with respect to Earth as a 

function of 6 ^ and 6 . 
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Fig. 4. 9 Rocket velocity magnitude with respect to Mars as a 

function of and 6. 



iJl- 



sin 0 



sin 0 



1-e cos 6 1 - e cos 0 , 



+ 



sin 



•1 / e - cos 0 
(l-e cos 0 ^ 



-sin 



/ e- cos 0 , 
Vl-e cos 0 



= U {-f n/GM (4-24) 

I Si 

where GM = 0. 01716 if a is given in AU. tv. is in days (sidereal). 

Values for a large number of trajectory time of flights v/ere 
calculated and the results plotted in Fig. 4.10. Fig. 4.2 illus- 
trates pertinent factors in the foregoing equations . 



4.12 Angular Separation of Earth and Mars 

Having found the time of flight for a given trajectory we may 

compute the angular .separation of Earth and Mars at the instant 

of firing (0 - 0 ) . Since Mars is in motion about the Sun(0 - 6 ) 

° o s os 

is a quantity which must be known in order to fire at the correct 
time to produce an intercept. It is evident from Definition Summary 1 
that 



+ ir". 



(4-25) 



/■.’7 



It may be shown ' ' that 



in 



(l-e cos 0 ^ 



9"' T ,i-e2)S/2 



(4-26) 



= constant (l-e cos 



(4-27) 



Since ^Q^^is a function of ©Q^ve are inexorably engaged with 

a quantity v/hich varies not only ivith real time but also with ang- 

0 

ular position. A. convenient method of expressing is fo plot 

00C^ as a function of time; then, knov/ing t^ and choosing t^ we can 
obtain the angle 

^ * tf — 6 - 0 A 

Oo f s A 






(4-28) 
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Fig. 4. 10 Time of flight as a function of 6 and 6 



and 



<«o 



«s> ' <«o - 



- ®A> 



(4-29) 



The aforementioned plot appears as Fig. 4. 11 . . To assist 
in the logical selection of t^ an overlay ( Pig. 4. 12; in en- 
velope attached to back cover) has been constructed. The basis 
of construction and description of use is given in Information 
Summary 4.1 



60 



MARS ANGLE IN DEGREES FROM EARTH APHELION 




F’ig. 4.11 Mars angular position measured from Earth aphelion 
as a function of time 
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V 



Fig. 4.12 (Overlay#!) 

Earth angular position measured from earth 
aphelion as a function of time 



(In pocket on back cover) 
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INFORMATION SUMMARY 4. 1 



Construction and Us e of Ove rlay #1. 

Assuming the "state of the art" of interplanetary travel is 

such that an Earth-Mars trip is planned to occur about 1980, 

2 54 o 

we note ' that on January 1, 1979 Earth is 177 ahead of her 

aphelion (see insert on Fig. 4.12 ). Mars is 2^ behind 

Earth's aphelion. In other words Mars is 181° ahead of the 

Earth. During the succeeding 1.135 sidereal years Earth 

overtakes Mars, and*they are in conjunction 227° past Earth 

aphelion (about February 17, 1980). Hence from January 1, 1979 

to February 17, 1980 ( 9^ ■* ®g) takes on all values from 0 to ir , 

With all pertinent parameters of both Earth and Mars orbits 
54 

known , Fig. 4, H and 4,12 were drawn so that information for 
the period in question would be easily accessible. The as- 
sumption has been implicitly made herein that Mars aphelion 
processes only slightly about the Sun (a few seconds of arc in 
about 10, 000 years). 

(a) Basis for Construction of Overlay 

The angular position of Earth is indicated for January 1,1955 
and will repeat every sidereal year. Mars angular position is 
marked on graph 4.11 for January 1, 1955, but its period is not 
an even multiple of years, hence its angular position for any 
given time must be found either from an Ephemeris or by the 
following method: 



(Page 1 of 3) 



Divide the difference in time between January 1, 1955 and 
the time for which Mars angular position is desired by the 
number 1. 88 (Mars period in sidereal years). Slide the mark 
for January 1, 1955 on the overlay to the right of the mark for 
January 1, 1955 on the underlying graph an amount equal in 
sidereal years to the remainder in the foregoing division. The 
real time scales have now been superimposed so that the 
angular positions of Earth and Mars are coincident on the time 
scale. The values of 6 ^ and 6^, rnay be picked off the curves 
and ordinate scale. If a time other than January 1, for a given 
year is desired^ make auxilliary pencil marks for Earth position 
and use that lather than January 1, as a reference. 

(b) To find ( 6 ^ - 0^) for a given t^ 

During the period of interest all values of ( 6 ^ - 6 ^) from 
0 to ;r occur, hence set the overlay as before and read t^ corre- 
sponding to the given 6 and 6 as a decimal fraction of a year 

o s 

after the chosen reference time. 

(c) To find t for a given 6 and 6 

o ° o s_ 

This problem is somewhat more complex. We are primar- 
ily concerned with a real time firing time in late 1979. This 

is true because it is about this time that (6 - 6 ) takes on 

o s 

values in which v/e are interested (note Section 4 J2 and Fig. 

4.11 and 4.12 ). Fortuitously; Mars' angular rate in late 

1979 is approximately constant at 0.47 degrees per sidereal 
day. If we multiply 0.47 times the t^ this determines a rough 
value of (6 - ®a)- Subtracting (6 - ^ from (6 - ^ a ) pro- 

vides (6^ - 0 ^) for which a t^ can be found as explained in part 
(c). By nothing the angular positions of Earth and Mars and by 
moving a straight-edge up or down the superimposed curves 
one may obtain a more accurate value of (6^- 0g) for the given 
data. 



(Page 2 of 3) 




(d) To find t for a given t„, 0 and 6. 

O I O 

January 1, 1955 was chosen as reference time because 
accurate values were available to the authors. 

The method described herein is a rough graphical method 
of finding required quantities for preliminary work. Actual 
calculation would of necessity be made by computer. 

(e) Notes 

(Page 3 of 3) 



4.13 V 0 ^ and ^ as functions of (9^ - 9^) 

Since we now may find (0 - 0 ) for any given trajectory, 

o s 

it appears pertinent to note any relations which may exist 
between this angle and the velocity required to fly the trajecto- 
ry. From an inspection of the basic equations in Section 4.8 
it is obvious that the relationship is quite complex. For this 
reason sind were plotted graphically for a large 

number of trajectories in a manner analogous to that of 



Preston-Thomas 
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The minimum ^ was found to be the 



envelope of the curves as 0^ took on various values. (See Fig. 

4.13 ) The minimum was found to be identically the 

curve for 0^ = 0°. (See Fig. 4.14 ) This is not inconsistent 
since the first case represented minimum velocities for a 
non-maximum radius vector of an ellipse; whereas the second 
case represented the minimum velocities that exist at aphelion 
(maximum radius) of an ellipse. In both cases the Hohmann 
Ellipse (0^ = 180°; 0^= 0°) provided the absolute minimum 
velocities . 



4.14 The Rocket Plane (A plane) Concept 

Let us define the plane containing Rq 0 and V 0 ^ at t = t^ 
as the rocket plane (A plane). If the rocket were not perturbed 
during its flight it could remain in this A plane (this being a 
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Fig. 4. 14 Rocket velocity magnitude with respect to Mars as a 

function of 6 -6 



44 

property of the motion of a body in a central force field ) . 

The A plane need not necessarily, and generally will not, lie in 
the Earth ecliptic plane since Mars passes through the Earth 
ecliptic plane only tv/ice every 1.88 sidereal years. Unfortu- 
nately, for calculation purposes, the orientation of Vq^ in 
Solar System space depends on the real time firing time and 
time of flight. Hence every possible trajectory requires a 
complete individual solution. The labor involved is so moun- 
tainous that solution by high-speed digital computer is manda- 
tory; and even these solutions must be based on a selected 
future firing time. 

If we can show that the displacement of the rocket normal 
to the A plane is negligibly small we may be justified in saying 
that for practical pusposes the rocket remains essentially in 
the rocket plane and our problem reduces to two dimensions. 
Neglecting the possible perturbations of bodies such as meteors, 
distant planets, asteroids, etc., the bodies causing displacement 
of the rocket from the A plane are Earth, Moon, and Mars . In 
an effort to show the minor nature of the aforementioned dis- 
placements calculations were made on a single "typical" 
trajectory. This trajectory was chosen because =10, 000 

ft /sec, which is a velocity probably within the capabilities 
achievable in the near future Since we are concerned with 
the period from late 1979 to early 1980 the firing time is es- 
tablished by ( 0 - 0 ) to be about December 6, 1979 Note 

Fig s . 4,15 and 4.18, 

Referring to Derivation Summary 4.1 the cumulative dis- 
placement of the rocket from the A plane was found by ac- 
cumulating the values determined from Equation 11 over the 
trajectory traveled by an unperturbed rocket(see Fig 4.19). 
Certain errors are inherent in the calculations but will be seen 
to be negligibly small. The arc length of the trajectory 
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ROCKET ELLIPSE DATE OF FIRING 




PERIHELION 



Fig. 4. 15 Angular quantities of a representative 
trajectory discussed in section 4.14 
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Fig. 4. 16 A typical trajectory showing the effect of tilt to account for 
non- CO- planarity of Earth and Mars motion 
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AN IMAGINARY CYLINDER 





Fig. 4.17 Assisting supplementary view to figure 4.16 



71 



< 




lO 



72 



Fig. 4. 18 Plan view of body motions for trajectory discussed in 

section 4.14 
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DERIVATION OF THE EQUATION FOR THE DETERMINATION 
OF THE DISPLACEMENT OF THE ROCKET NORMAL TO THE 

"ROCKET PLANE" 



Four relations are available from consideration of the geo- 
metrical properties of the trajectories or orbits. Two assumptions 
are made which greatly simplify the solution. The error introduced 
by these assumptions is on the "high" side and will be shown to be 
negligible (Fig. 4.16 refers). 

'a ' <‘max’A ^ 

(A + - 6(5 > 

z . = R^ . sin i . 

A CA A 

Z(^= sin 

From (3) and (4) 

. -I 

1 . = sin 

A 

(a) Basic Equations 

To aid in clarification of the development that follows and 
the assumptions made, the following list of symbols is presented: 



<1 



R 



CA 



sin 



(5) 



( 2 ) 

(3) 

(4) 



^^orth^j 



Magnitude of that component of normal to 
the "rocket plane". 



DERIVATION SUMMARY 4.1 
(Page 1 of 4) 




i. Inclination of indicated quantity with reference 

J 

plane, or with ecliptic plane if no reference 
is given-positive direction in the conventional 
sense . 



(i ). Maximum inclination as defined above, 
max j 



z. 

J 



A 



£ 



Displacement of indicated quantity normal to 
indicated reference or to ecliptic plane if no 
reference is given. 

Angle as shown in Fig. 4.16 

Eccentric anomaly as shown in Fig. '±. 16 

Increment (Magnitude as defined). 



(b) Definitions 

Referring to Fig. 4. 16 and4.17the two assumptions made are 
( 1 ) that the "rocket ellipse" is almost circular about its 
center (C). 

( 2 ) that the line CC in Fig. 4.16 is negligibly small. 

(c) Assumptions 

To derive the required equation we solve (3 ), then 

( 4 ) for the final condition 6^. We then find i^ from 

( 5 ) and (i from (1 ). Having discovered (i ) . 

for a given trajectory we are now able to write i^ as a function 
of (Equation 1 ). 

Any [ z ] may be found thusly: 

DERIVATION SUMMARY 4.1 
( Page 2 of 4 ) 



From Fig. 4.17 ; 



and 



obviously 



f ^Aplane^ ^ ^C-A ^A^ 3-^(|5 



Z(j = [r<5 sin 



^ Aplane ^ ^Aplane^ at^ ^CJ 



now 



[ ‘di 



-1 



Apian e 

measured from 
rocket 



sm 



(Ml 



Aplane 



'd- ■'a I 



then 



and 






ISorthld ° '■a> 



[Zdl 



Aplane 



I'd' ’‘a 



1 *^a 1 Aplane Wue to Mars) 

= Lf^Ithd ( 5„2 

if 6 t is in sidereal days 

8orth *" ft/sec^ 

6 Z^is in miles 

DERIVATION SUMMARY 4.1 
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( 6 ) 



(7) 



(8) 



(9) 



(10) 
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A similar derivation may be made for the cis -lunar 
Barycenter . 

(d) Derivation 
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traveled by the rocket is approximately 354, 000, 000 miles; the 
cumulative displacement from the A plane is approximately 5, 500 
miles: a ratio of about 1. 55 to 100, 000, certainly negligible. 
Further, the calculations show that the displacement due to the 
mass attraction of Mars is only 33 miles, the remainder being 
the effect of the Earth-Moon mass attraction. About 3,250 
miles displacement occurs in the first five days of travel; hence, 
if a computer program were established to counteract this, the 
displacement could be reduced even further. 

The obvious advantages of reducing the problem to two 
dimensions are many -fold. These include (a) reduction in 
weight and complexity of the guidance mechanism; (b) simpli- 
fication of component construction; (c) savings in fuel; (d) im- 
provement of payload ratio; (e) saving in material, labor, and 
cost; (f) contraction of production time-tables; and many others. 



CHAPTER 5 



DEPARTURE AND ARRIVAL 
ORBITS 



5.1 Introduction 

The capability of a rocket to achieve required velocities de- 
pends so critically upon the energy carried within itself that 
every effort must be made to minimize these energy requirements. 
This section investigates one aspect of this minimization -- the 
selection of departure and arrival orbits about the respective 
planets to minimize kinetic energy, hence velocity, that the 
rocket must provide itself in order to fly a given trajectory. 

Fig. 5.1 illustrates the relationships of pertinent velocities. 
For complete ctefinitions see Definition Summary 1 




79 



5. 2 Optimum Orbit Radii 



In order to fly any given trajectory, or ellipse, a velocity of 
the rocket with respect to the Sun (Vq must somehow be attained. 
Referring to Section 4.9, we note that this is the velocity vector 
required for an intercept of the orbit of Mars,^ treating both Earth 
and Mars initially as massless points. However, this is obviously 
not the case. Referring to Fig. 5. 2 we can say, figuratively, 
that we must climb up the side of the Earth potential "well", over 
the intervening "hump", then down into Mars potential "well. 
Jupiter's potential "well" is shown for comparison. From Fig. 5.1 
we note that the rocket already possesses a certain velocity with 
respect to the Sun by virtue of the Earth's orbital motion, thus the 
rocket must be provided the velocity, V 0 ^. In order to "over- 
come" the effect of the earth gravity field, the rocket must be 



provided a greater velocity, V 



(0A)g- 



The rocket also possesses 



a velocity, V 0 by virtue of being in orbit about the earth, 
therefore it remains for the rocket to provide by its own efforts 
an increment of velocity Initally Vq ^ will be along the 

dotted line labeled in Fig. 5.1, however, the earth gravity effect 
decreases rapidly so that ^ lies almost identically on the solid 
line. In addition the ratio Irm /R/ 0 /ii| is of the order 10 , hence 



the displacement of the rocket from the earth is negligible 
for purposes of calculating initial velocity requirements. 



Regardless of validity, if Vq ^ is known we may obtain 

relations between V. and r/jN , This has been done in Derivation 

me s 

Summary 5.1. Although the relations have be^n developed for the 
planet Earth the derivation is perfectly general and applies to any 
body; with different constants. 



Parts (b) and (c) of Derivation Summary 5.1 develop and 
prove the optimal relationship of and 1 * 0 ^. Part (d) expresses 

the relation non- dimensionally . 

Pertinent ranges of vs T 0 ^ have been plotted in Fig. 

5. 3 and the non-dimensional function is shown over an applicable 
range in Fig. 5.4. 
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Fig. 5. 2 Velocity of escape" from the Sun and some planets 

as a function of distance 
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Fig. 5.3 Optimal relationships between and for 

Earth and Mars 
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is the frequency of occurence of the rather special Earth-Mars 
configurations required bythe Hohmann ellipse. This, among 
other considerations has been taken up in preceding sections. 



5, 3 Use of Optimum Orbit Concept 



Almost paradoxially., the optimum orbit concept is of 
particular use in the selection of non-optimum orbits or tra- 
jectories. As an example, assume that a rocket has been 
developed with a capacity to provide itself with a of 

approximately 10, 000 ft /sec The optimal orbit radius corre- 
sponding tc 10,000 ft/sec, as an optimal velocity, from graph 5 3 
is approximately 27, 500 miles above Earth center. If 27„ 500 
miles is unacceptable due to other considerations we can reduce 
(from graph 5,4 ) this radius by 42 . 5% simply by allowing a 2% 
increase in the required velocity. In other words, for 
= 10, 200 ft /sec r = 15, 800 miles, or approximately 12, 000 
miles above the Earthy's surface. It must be borne in mind, 
however, that these relations exist only for that trajectory for 

which (V. ) , = 10,000 ft /sec 

me opt 

Applying the foregoing reasoning to the "Hohmann" type 
ellipse -(V. ) , = 6, 818 ft/sec; (r * ) , = 57, 300 miles - we 

may reduce r^^ to 33, 000 miles by increasing to 6, 936 

ft /sec. Hence, if we can provide the energy to increase 
by 118 ft/sec we not only gain the obvious advantages of a lesser 
r but also the flexibility inherent in a greater freedom of 
choice in the angular positions of Earth and Mars. This has been 
discussed more fully in Chapter 4 

5 4 The Least Orbit Concept 

If the problem is posed that = 7, 500 ft /sec. , what is 

the least orbit the rocket may depart from and still travel a 

trajectory to arrive at Mars ? This obviously is not an optimum 

orbit problem as previously posed. For > ^00 ft /sec 

®s^opt " miles (graph 5.3). For = 6, 818 

ft /sec (the minimum "Hohmann" ellipse velocity) the ratio 

/(V. ) , = 11 and from graph 5.4 the ratio r^ , 

me opt & ©s ®s opt 

= 0.327 and r^ = 18, 750 miles The limiting restriction herein 
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DEPARTURE AND ARRIVAL 
ORBITS 



Optimum Planetary Orbits 



In order to "escape” any bodies mass attraction a vehicle 
must attain kenetic energy at least equal to the potential energy 
inherent in the field or fields to be "overcome"; hence 



KE = m 



or 





ilV 


dR 


(lA) 


00 


GM. 








3 

o 


dR 


(IB) 


■ R 


R 







KE 

Following the convention Q = |q| we integrate 
KE 



and 



m 



GM. 
] 

^ A 

JA 



( 2 ) 



but 



KE 
m * 



1/2 V 



je 



(3) 



hence 



V. 

je 



2GM. 

V 

63 



(4) 



It has been shown that a rocket in circular orbit about a 
body possesses a so-called circular velocity. 



DERIVATION SUMMARY 5.1 
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'GM. 



V. = 

]s 



R. 



hence 



jA 



V. = V2V. 
je js 

From Fig. in the case of the earth and a rocket 
(®A)g me ©s 



(5) 



( 6 ) 



(7) 



Further, by equating the kinetic energy requirements for 
the rocket we have. 

1/2 m.V .^+l/2m.V ^ 

^0 ^ ^0 ^ A0 (0A)g 

( 8 ) 

where the first term represents the kinetic energy requirement 
to achieve the velocity fhe second term preresents the 

kinetic energy requirement to "escape" the earth gravity field; 
and the third term describes the sum of the first and second. 



From (8) we have 






Hence, substituting in (7) and rearranging 
V 



me 



= nA? 



0 A ^ ^^0 s 



V, 



0s 



(9) 



(10) 



and if initially R . = r then 

A s 



V 



me 






‘'0 A 



*^0s 



"0s 



(11) 



(a) Development of V in terms of rocket orbital radius 
^ me 
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In order that r^^ be the optimum radius on a kinetic energy 

basis V. must be a minimum for r/rv =(r/i, ) , which is the 

me tus CPS opt 

equivalent of the conditions »■ 

6V. 



( 1 ) 



me 



6n 



= 0 



0s 



6^V. 

( 2 ) — > 0 



( 12 ) 



(13) 



6r, 



■0 



To satisfy (1) we take the first partial derivative of (11) 
and equate to zero. 

GM^ 



GM 



0 



6V. 

me 






2 

"“es 



--;2GMe 

V/TN A » 



+ 1/2 






^0 A n 



0s 



“■®s 
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To satisfy condition (2) above we take the second partial 
derivative of with respect to r^^ and evaluate at r^^ 

= (r 0 g) ^ to show it is indeed positive. 
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(c) Proof of Minimality of (r^^)^^^ 

Dividing (11) by (17B) and noting (17 A) we have 
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(d) Non-dimensional Relations 
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CHAPTER 6 



THRUST AND ENERGY CONSIDERATIONS 
A, THRUST REQUIREMENTS 

6. 1 Introduction 

It IS of interest to examine the thrust requirements and 
payload ratios for a rocket departing circular orbit about earth. 
Such figures appear to be more or less well known for the 
journey from earth's surface into orbit (as for earth satellites). 
However, the orbit-to -orbit requirements are not commonly 
available, although easily calculable. This information is con- 
sistent with the optimal orbit concept developed previously and 
verifies the validity of such an approach for the conditions 
assumed. 

Although the data have been based upon the departure and 
arrival orbit concept, an extension to surface departure can 
easily be made by extrapolation to lower altitudes , The penalty 
of departing a lower altitude will be evident after an investigation 
of this section 

6. 2 The Optima] Payload Ratio 

The change in velocity necessary to affect transfer from 
orbital motion about earth to departure path has been designated, 

. This change in velocity is accomplished by applying thrust 
for a controlled period of time. If the power-on period is suf- 
ficiently brief the maneuver trajectory can be neglected and the 
transition from circular orbit to departure path thought to occur 
instantaneously at tangency. The justification of this assumption 
will be discussed shortly after the development of the equations. 

The commonly used expressions for an idealized single 
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step thrust study of this nature are derived from the equations of 

6 3 

motion written along the orbit flight path. 

m-^ = T - D - Mg sin 6 

In circular orbit, however, for the above assumptions, D 
is presumed negligible and 0=0. 



, dv 
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M t , 

e_ for 0<t<t 

t^ dt -^b 



where = initial gross mass 
Mp = mass of propellant 
t^ = time of burning 



A constant thrust device is assumed of specific impulse, I. 
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Integrating between t = 0 and t = 



V - V = V. = - Ig In (1 - 6)= - Ig In (PR) 
t, o me . & ' , & ' 

b 



PR = Payload Ratio 



M - M 
o P 

M 



It is to be remembered that this, V. , applies only to 

mc0 '' 

the departure maneuver, as does PR0 • For entry into circular 
orbit about Mars, a similar expression is obtained. The total 
payload ratio ( PR.|.) for an orbit -to -orbit maneuver, ignoring 
corrective thrust requirements , is 



PRt - PR® PR cT 



These equations have been expressed in tabular and 
graphical form for various values of specific impulse^ which 
has been chosen as the propellant parameter. The subscripts 
in the headings of Table 6.1 refer to the three assumed values 
of specific impulse. Table 6, 1 for example, lists Payload ratios 
for various sample trajectories developed by the methods of 
chapter 4 . The importance of improved propellant performance 

is noted. Although the corrective thrust requirements have not 
been included, reference 47 
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TABLE 6.1 

SOME REPRESENTATIVE TRAJECTORIES WITH CORRESPONDING 

PAYLOAD REQUIREMENTS 
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has indicated that a suitable payload ratio for corrective thrust 
requirements is .79. If the above figures are altered by this 
multiplier, it is seen that the most optimistic Payload Ratio to 
be hoped for is .11. The calculations of Table 6.1 presume de- 
parture and arrival orbits at "optimal" altitudfe as defined 
earlier. Since these altitudes inevitably vary for various 
trajectories, little correlation is to be expected between Payload 
Ratio (PRrp) and angular separation at instant of fire (9^ - 6^) 
other than to indicate a lower bound for certain specified para- 
meters , 

6. 3 An Interpretation of "Optimum" Altitude 

At first glance, it appears that the higher the optimal 
departure and arrival orbit altitudes, the more favorable the 
total Payload Ratio. To invalidate this idea as well as to demon- 
strate the advisability of employing optimal altitude two tra- 
jectories have been selected for study: one is the familiar Hohmann 
semi-ellipse trajectory: the other merely a representative one 
selected for contrast. Both are defined by the * in Table I. 

If the PR^ and PR(J* for these two trajectories are 
calculated as functions of altitude (including optimal altitude) on 
a non-dimensional basis, as in Fig. 6. land 6. 2 it can be seen 
that a peak payload ratio is achieved at precisely (r )qp'p> A slight 
fall -off occurs at values of (r f > 1. Thus, the choice of 

S Qprj. 

optimum altitudes, in the sense defined, insures a maximum 

total payload ratio for any given specification of V. ^ and 

inc w 

^incCf ■ great saving in Payload Ratio is noted when one 

departs from circular orbit altitude for ^ .5. 

<''s'0PT 

6. 4 The Assumption of Impulsive Thrust 

A justification of the assumption of impulsive thrust is 
dependent upon a comparison of burning time and period of orbit. 
Obviously, for a low thrust device of considerable duration, the 
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Fig. 6.1 Payload ratios VS. orbit altitude for the Hohmann ellipse 
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Fig. 6.2 Payload ratios vs. orbit altitude for a representative trajectory 
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PERIOD-MINUTES 




Fig, 6.3 Orbit pt^riod as a function of altitude above planet center 
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assumption is not valid. 

The period of orbit for both Mars and Earth is plotted in 

Fig. 6.3 as a function of distance from Earth and Mars center. 

If we assume that instrumentation is capable of withstanding an 

axial load of lOg upon departure from Earth o^bit the minimunci 

0 3 

acceptable burning time can be determined . Careful design 
can then be employed to produce a favorable burning time. The 
minimum burning times as functions of orbital period for optimum 
Earth and Mars orbit are listed in Table 6. 2 for the trajectories 
previously considered: A specific impulse of 350 seconds is 
assumed. 
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TABLE 6.2 

A COMPARISON OF ORBITAL PERIOD AND BURNING TIME 
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B. ENERGY CONSIDERATIONS 



6 . 5 Introduction 

It would seem appropriate at this stage to mention opti- 
mization of the trajectory. The term, optimization, however, 
invariably leads to ambiguity and requires explanation. Do we 
desire optimization with regard to time, fuel expenditure, or 
initial velocity?-- or possibly cost? Previously, an attempt was 
made to minimize the initial velocity requirements, but it is 
clear that this does not necessarily represent the optimum as 
far as other variables are concerned. 



An overriding consideration during early interplanetary 
flights will undoubtedly be minimization of fuel expenditure. By 
minimization of fuel expenditure, we imply a determination of 
not only the minimum duration of powered flight but also the 
direction in which this power or thrust is applied. Accordingly, 
these are all considerations in the minimization problem. It is 
the intent of this section to describe the problem in general terms 

and summarize, without proof, the findings of other authors who 

^ X 1 XU XU- -1 u- X 34, 35, 36, 37 

have written at length on this particular subject. 



6.6 General Theory 

Let us first consider the problem of optimum transfer (with 
respect to fuel expenditure) between two terminals lying in a 
plane. If we consider an inertial orthogonal coordinate system, 
the equations of motion for the rocket are 



t • 

X + 



* • 
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G - 
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-cosa^ 


dM 
d t 


(6-1) 


G - 
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dM 

dt 


(6-2) 
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(6-3) 



Z + G 

z 



-cos 7 



c 

M 



dM 
d t 



where G ^ G , and G are the components of the resultant at- 

tractive force per unit mass and a, ^ and y are the angles 

between rocket thrust direction and coordinate axes. The ex- 
c dM 

pression is the well known simplified term for rocket 

thrust and is defined in the section describing thrust consider- 
ations . c is the effective exhaust velocity and M is the mass of 
the rocket at any time t. 

G , G , G depend explicitly on x, y, z, and t only, 

X y z 

Squaring and adding 6-1, 6-2, and 6-3. 



9 *♦ o 99 9 n HM 9 9 9 9 

[ X + G^] Y+Gy] + [ Z+G J = [ ^ ^] [ cos^ a+ cos^^+ cos^ 7 ] 

(6-4) 

1/2 



Integrating, 



M 



c log 






n 99 p 99 p 

[X+G^] +[Y +Gy] +[Z + GJ^ 



dt 



It is assumed that the boundary or terminal conditions are 
specified 

X^ (t^) = X^ 



) = X 


X (t ) = 
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Xj(tj) = X. 
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Y (t ) = 


Y 
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Z (t^) = 


Z 

0 


Z^(t^) = z 



Yi (t^) - Y^ 
(tj) = Z^ 



The objective is to select the functions x(t), y(t), and z(t) so 
that the fuel ratio is a minimum (i, e., we desire initial fuel 
mass, as small ^ as possible and still attain our objective 
subject to the boundary conditions). 

The minimization of the above integral is a problem in the 
calculus of variations. Its solution is highly detailed and is not 
presented here. The conclusions and interpretations of existing 
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solutions are helpful and interesting, however, and are stated 

37 

here without proof 



These conclusions are due to D. F. Lawden, one of the 
more active writers in the field. 

Let us consider a rocket orbiting the earth circularly. The 
optimal path is desired for travel from this orbit to entry into 
circular orbit about another planet. The conclusions of Lawden 
indicate that an optimal path can only be followed in a vacuum by 
the application of certain impulsive thrusts from the rocket motor, 
the rocket coasting freely under gravity during the intervals 
between such thrusts. These coast periods are essentially null- 
thrust arcs and are mihimal in the sense that neighboring arcs 
can be followed only by an actual expenditure of fuel. To establish 
an absolute minimal mode of transfer for a definite planet con- 
figuration, cases for various numbers of impulses must be com- 
35 

puted and the minimum selected. Lawden suggests that two 
impulses usually lead to the optimal solution. It should be re- 
membered, however, that there will be a large number of optimal 
solutions corresponding to different planet positions. Only one of 
these will be an absolute optimal solution for all planet con- 
figurations. The absolute optimal solution for an Earth-Mars 

flight has been known theoretically for some time if certain 

22 

fundamental assumptions are made. The analysis is that of 
Hohmann, who concluded in 1925 that the problem of transferring 
a rocket from one circular orbit into another coplanar with the 
first and about the same attracting body would be most economical 
if effected along an elliptical trajectory tangential at its apses to 
both circular orbits. Fuel would be expended rapidly at each apse 
to carry the rocket into and out of the transfer ellipse. A state of 
free flight would exist between the two impulse trajectory. This 
classic development, the so-called Hohmann ellipse, has been 
verified by Lawden and is consistent with the remarks above. Note 
that the optimum path as developed and computed in this paper on 
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the basis of a minimum initial velocity is in fact the Hohmann 
ellipse. The importance, however, of knowing the optimum for 
other positions of Mars and Earth so that flight can be initiated 
often, is not to be underestimated. 







Decelerating 

Impulse 



Accelerating 

Impulse 



Fig. 6.4 The Theoretical Mode of Transfer 



6. 7 Rocket Thrust Direction 

Ordinarily, the mode of directing the thrust which results 

in maximum total energy is found to be alignment of the effective 

exhaust velocity with the velocity vector of the rocket (i.e. align 

the direction of thrust with the tangent to the trajectory. This 

will ensure that the rate of increase of total energy is maintained 

at its maximum value at every instant. This program will be an 

economical one since the object of the maneuver is to raise the 

total energy as rapidly as possible in order to achieve escape 

34 

and sufficient velocity. However, it has been proved that this 
program may not be an absolute optimal method, although close 
to it. Under some circumstances, when departing a circular orbit, 
it is not economical to maintain the rate of increase of energy as 
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its highest value over the entire trajectory of escape. The 

author claims that after having escaped, it is advantageous first 

to direct the thrust so as to oppose the motion. The rocket then 

drifts toward the center of attraction along an elliptical arc. At 

its closest point of approach, an impulsive thrust is generated in 

the direction of motion. The velocity is now large, and the total 

energy increases at a rapid rate. The rocket is transferred into 

a hyperbolic orbit along which it coasts to its destination and the 

overall fuel consumption is found to be less than if escape had 

been achieved by application of a single tangential thrust. The 

3 8 

reader is referred to reference for information as to how 
advantage can be taken of such maneuvers for general conservation 
of energy. A good approximation to the optimal programming of 
thrust direction under departure from circular orbit may be 
achieved by maintaining alignment between directions of thrust 
and motion, 

51 

Other writers consider it an empirical rule that a journey 
from one to another stable orbit round a central body will be 
completed in the shortest time only if power is applied in two 
impulses, one at beginning and one at end of journey. 
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CHAPTER 7 



NAVIGATION IN SPACE 

A, POSITION IN SPACE 

7. 1 General Discussion of Problem 

Mid-course interplanetary navigation differs quite markedly 
from terrestrial position determination due primarily to the need 
to obtain a fix within a volume instead of upon the surface of a 
sphere (Earth) of known radius. Terrestrial positions can be 
obtained by simple spherical trigonometric relations and can be 

9 

further simplified to a problem in plane geometry and time. 

Although it IS obvious that the problem of Position in Space 
is three dimensional^ several writers have postulated that a suf- 
ficiently accurate solution can be obtained by considering the 

planets to be coplanar and working with only two dimensions in 

3 68 

the plane of the ecliptic. ^ Part of the argument for this 

approach is found in star tracker limitations, errors in astro- 
nomical measurements, the assumption that terminal guidance 
will correct for errors incurred during mid-course transit, and 
the relatively small angles of inclinations of the planet orbits 
from the ecliptic plane, 

7. 2 Factors Indicating Necessity for "On Board" Automatic Solution 

If we postulate that accurate positional data is needed during 
mid-course transit, several factors indicated that the problem 
must be solved on board the vehicle and by automatic computers 
of some sort. 

1, Due to high vehicle velocities (in the neighborhood of 
1, 200 miles per min) time delays of any sort should be 
avoided wherever possible. This eliminates manual 
computation and sight taking in the vehicle itself or on 
the earth. 
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2. Radio frequency transmission in space has a rather 

0 

sizeable time delay, (A 10 mile round trip signal requires 
approximately 18 minutes) 

3. Radio frequency signal interference in space is un- 
certain and apparently unpredictable Vhen sent to or 
from Earth due to the Earth's atmosphere and extraneous 
radio frequency waves in space itself.^ 

4. Power requirements for dependable radio trans- 
missions during mid-course appear to be far too great for 
installation on any vehicle within the near future , 

5. Optical tracking of the vehicle from earth is not possi- 

69 

ble with any foreseeable astro telescopic equipment. 

6. Radar tracking of the vehicle from earth is ruled out 
by 2 and 3 , 

7. 3 Selection of a Reference Coordinate System 

Celestial bodies may be located with respect to several 
different reference systems which are shown in Fig. 7.1 and 
Table 7.1 . For mid-course interplanetary position de- 

termination two of these systems can be discarded immediately. 
These are the Geocentric Horizon and Equatorial. Coordinate 
systems using Galactic quantities are not particularly convenient 
for travel limited to our solar system. Of the remaining only the 
heliocentric ecliptic system will be considered for use in this 
paper. This system has the following features which lead to 
this choice. 

1. The reference zero latitude plane includes the Earth 
and Sun at all times. 

2. All other planet's orbits lie close to this plane, 

3. NumerousEphemerides arc available based in this 
system simplifying programming of any memory drum for 
a computer 
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Fig. 7.1 Reference coordinate systems 



TABLE 7.1 CCX)RDINATE SYSTEMS 
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7 . 4 Establishing the Reference Coordinate System 



It appears that the reference system orientation will have to 

be obtained by use of the fixed stars, an extremely sensitive gyro 

71 

inertial platform, or a combination of the two. The use of 
gyroscopic elements by themselves for establishing and main- 
taining a chosen reference coordinate system to a sufficient degree 
of accuracy to determine accurate angular measurements and rates 
in space appears to be beyond engineering capabilities at present. 
All instrument gyros have drift uncertainties to some extent and 
very definite lower limits of sensitivities. To be useful for ac- 
curate interplanetary guidance the threshold of sensitivity would 
have to be in the order of lO"* g. Since it is necessary to con- 
template transit times of over 100 days, it appears that star 
trackers locked on the fixed stars will have to be used to a cer- 
tain degree in either assisting a gyro system or replacing it 
entirely. At certain periods of the transit it may be advantageous 
to use gyroscopic elements as a memory circuit or coast device 
while a star is temporarily obscured or its line of sight becomes 
nearly coincident to that of a stronger source of light such as a 
planet or the sun. It also appears that during periods of radical 
maneuvering or high thrust a gyroscopic device may be needed 
to insure holding track on a desired star. 

If we postulate the use of photoelectric star trackers, es- 
tablishing the desired reference coordinate system can be ac- 
complished by using only two fixed stars approximately 90° apart 
in the same meridian. However, three stars near quadrature to 
each other, or three stars near the chosen zero latitude spaced 
approximately 120° apart, will establish a more stable and secure 
reference. Monitoring three fixed stars would allow for the 
possibility of one of them becoming momentarily obscured or 
temporarily loosing tracking signal. Section 8. 64 discusses 
the accuracy requirements for such systems. 
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7 5 Determining Celestial Latitude and Longitude of the Vehicle 



After establishing the reference coordinate plane and meridian, 
the celestial latitude and longitude of the vehicle can be determined 
from the line of sight to the sun as seen from the vehicle. This 
is illustrated in Fig. 7,2 and as shown, the veliicle heliocentric 
latitude is equal to the negative of the angle between the reference 
zero latitude plane and the vehicle sun line measured in a meridian 
of the reference coordinate system. The heliocentric longitude of 
the vehicle is equal to the angle between the reference zero longitude 
and the meridian containing the vehicle sun line measured in the 
direction of rotation of the planets, plus 180°. 

7. 6 Determining the Radius Vector to the Sun 

After establishing the basic reference, and determining the 
vehicle latitude and longitude with respect to the sun, it is then 
necessary to obtain the distance of the vehicle from the sun. This 
can be accomplished using only one planet tracked with respect to 
the established reference system by a simple method proposed by 
Vertregt. 

We can obtain from a memory drum or an ephemeris of 
heliocentric ecliptic positions of the selected planet as a function 
of time; 

1 Longitude of the planet ( X ) 

2 . Latitude of the planet ( L^ ) 

3, Distance of planet from sun ( R ) 

To obtain the radial distance of the vehicle from the sun 
referring to Fig. 7. 3, we can proceed as follows: 



liO 



REFERENCE 





( 7 - 2 ^ 

Fig. 7, 2 Determination of heliocentric latitude and longitude 
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Fig. 7.3 Geometrical relationships determining radial 
dis+ance from Mie .sun 
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R 



(E-0) 



R 



Rp cos Lp sin [Ap - A(T.p)J 

' “= ^E-0) • - ^(T-p) >1 



(7-3) 



If the earth is one of the planets used for determining positon, 
and we use the plane of the ecliptic as the zero latitude reference 
plane, the solution is somewhat simplified due to the fact that 
earth heliocentric ecliptic latitude is at all times zero. (Fig. 7. 4 ) 
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Fig. 7.4 Geometrical relationships determining radial distance 
from the sun with earth as the reference planet 
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The determination of range of the vehicle from the sun can 
be made using only one planet after and are known. 

However, since we know appreciable errors exist in the values 
given for planetary positions as shown in Section 8. 62, the use 
of two or more planets should tend to narrow the error down. 

7. 7 Some Additional Methods for Position Determination 

A second method for angular determination of the sun and 

6 7 

planets is photographic v-.' .) • By the taking of simultaneous 

pictures of a given small area (such as 5°) centered about the 
selected body, the position of the body with respect to the 
background of recognizable fixed stars can be obtained by taking 
measurements from the developed plate. This does not appear 
practical even for a manned vehicle, due to time lag, the need 
for continuous or very frequent data to effect smooth system 
performance and, in the case of a manned vehicle, an extra 
burden upon the crew. 

Similarly it has been suggested that the angular measurement 
of the diameter of solar bodies might be used for range determi- 
nation. This could be photoelectric or photographic in nature. 

Both these systems do not appear practical except possibly for . 
terminal guidance where the diameter of the body being approached 
becomes appreciable. 

Further postulated ranging methods are (1) measurement of 

49 

incident solar radiation intensity, (2) measurement of reflected 
solar radiation from planets, and (3) radar ranging. 

Some basic arguments against use of the above methods for 
mid-course ranging are as follows; 

1. The variation in intensity of solar radiation does not 
seem to be great enough for any but the most coarse range 
determination. Using Fig. 2. 2 we see that for a trip 
between Earth and Mars this variation in radiation is 
approximately 1 erg per square inch per second for a 
change in range of one mile . 
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2 . Reflected solar radiation also has the same general 
gradiant as above in (1) plus the uncertain and possible 
widely varing amounts that might be reflected from a 
body do to changes in planetary atmospheric condition 
and orientation of lighted surface relative to the vehicle. 

3. Radar ranging has the same disadvantage as listed in 
Section 7. 2 

7. 8 Determination of the Shape of Actual Transfer Ellipse 

If we assume that after obtaining a series of acceptable 
vehicle positions we want to know the elements and orientation of 
the ellipse being traversed^ we can proceed after first converting 
the differences in latitude and longitude of three selected positions 
into angular displacements along their great circle locus ($2 - 0 ^). 
This can be accomplished by standard coordinate transformation. 

Determination of three of the basic elements of the unpowered 
orbit, 6 , e, and a, can be accomplished by the use of three se- 
lected positions of the vehicle with respect to the sun. For any 
reasonable degree of accuracy the time interval between selected 
positions may have to be several days , A method of mathematical 

reduction which can be easily programmed for a computer has been 

6 7 

suggested by Vertregt and is adapted herein to the nomenclature 
of this paper in Derivation Summary 7.1, 

Two assumptions are made in this method which are common 
to many proposed non-powered interplanetary mid-course tra- 
jectories . These are; 

1. Vehicle is following an unperturbed Keplerian ellipse, 

2. Sun is only gravitational force considered, and is one 
of the foci of this ellipse. 

From Derivation Summary 7.1 we obtain 



DETERMINATION OF THREE BASIC ELEMENTS OF ACTUAL 
ELLIPTICAL TRAJECTORY FROM OBSERVED DATA 



The following method of determining 6, e, and a is adapted 
from reference 67 to the nomenclature of this paper. 

Using a basic equation of any ellipse and referring to Fig. 7. 5 
we obtain the following: 



a (1 - e^) 
1+ e cos 6^ 



(5) 



a q - e^) 
^A 2 1 + e cos ©2 



( 6 ) 



a ( 1 - e^) 




(7) 



Combining (5) and (6) and (5) and (7) we get 




e cos 0, = R. + R . e cos B 
'■ ■^2 -^2 



2 



( 8 ) 




3 



(9) 



Rearranging (8) and (9) 




R 



A 



R 



A 



2 



1 



( 10 ) 



e 



R. cos 0, - R, cos 0 

-^1 1 -^2 



DERIVATION SUMMARY 7.1 
(Page 1 of 3) 
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Fig. 7.5 The rocket ellipse in space 
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From ( 11 ) we obtain 



Ra^(Ra 



I I 

R. ) cos 0, - R. (R. - R. ) cos 0„ 
2 ^2 ^ ^2 "^1 ^3 2 

+ R . (R. - R . ) cos 0„ = 0 



^2 ^^2 ~ 



0 „ = 0 , -t ( 0 „ - 0 , ) 



(12) 

(13) 

(14) 



II II 

(02 - 0 j^ ) and (^3 “ ^j^) have been determined along with R^ , 

Ra > by establishing the three instantaneous positions of the 

veliicle. ^e can now substitute (13) and (14) into ( 12 ) 



R. (R. - R. ) cos 0, - R. (R. - R. ) cos 

Aj Aj A3 1 Aj Aj A3 



V<«2- V 



. + R - (R . - R . ) cos 

^3 ^1 ^2 



0 ;.( 0 ;- 0 ;) 



= 0 



(15) 



Expanding the summed angle terms in (15) leads to 
R. (R. -R. ) cos 0, “ R. (R. -R. ) 

^2 ^3 ^ ^2 ^1 ^3 



I I I 

cos 0 ^ cos ( 02 “ 0 j)- 



I I I 

sin 0 ^ sin (02 - 0 ^) 

I I r 

sin 0 j^ sin (^3 - 0 j^ ) 



+ R» (Ra “ Ra ) 

^ ^2 

= 0 



r 

cos 0 j^ 



COS (02" 0J^) 
(16) 
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Solving (16) for 6^^ we get 



e; . tan-1 






I- 11 Ag) < ®2 - ®l‘>+ llAg'llA^- «A2> <®3 ' V 

(17) 





From (11), substituting the value obtained in (17) we obtain 





Ra - Ra 

® “ 1 t 

Ra cos 6, - R. cos 0„ 

Al 1 ‘is 3 


(18) 


Substituting the value obtained in (18) into (7) and rearranging 
leads to 




R^ (1+ e cos ©g) 
1 “ e 


(19) 
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«A3> <®2 ‘ ''a3<«Aj- «A2> ' V 




^ «A3> Aj- «A2> "i" <®3 * 





(7-5) 



(7-6) 



(7.-)7) 



7. 9 Determination of Orientation of Actual Transfer Ellipse 

Two additional elements of the transfer ellipse, and i, are 

required for determining the orientation of the vehicle trajectory 

with respect to the reference plane. These can be obtained by 

6 7 

applying the method suggested in reference using two positions 
of the vehicle. A rather obvious observation can be made that the 
further these two positions are from each other the more accurate 
will be the values of Q and i obtained. However, since corrective 
action is more economical if performed early in the flight, it is 
apparent that $7 and i should be determined as early as possible. 

If the trajectory is satisfactory at this time more accurate values 
can be obtained at the next determination. 




From the first vehicle position, as illustrated in Fig. 7. 6 we 
can obtain the following: 



is to NN 



A, 
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^1^0 ”^A sin 
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Fig. 7. 6 Orientation of transfer ellipse 



therefore 



sin L 



A. 



sin 1 = 



sin T7 



(7=8) 



1 



For the second position we can likewise determine that 



sin L 



sin 1 



sm[ T7^+ (02 ■■ 



Combining the above 



sin L 



sin r) 



sin L 



sin L 



^ = . ^ 

1 sin[r]j^+ (02’-” ^j)] sin 77, cos ( 0^ “ cos t]|^ sin(02“ 0j^) 



Multiplying by sin 17^ and regrouping leads to 



- tan 



sin sin (02 = 0j) 

_sin = sin cos (02 " 

2 1 



[7=9) 



can now be converted to the plane of the ecliptic. If we then 



subtract it from X we have fi. 



After determining we can then solve for the inclination of 
the orbit of the x'ehicle using equation (7-=8) 




(7-10) 




■J 



B. IN-FLIGHT CORRECTIONS 



7,10 Introduction 

It is considered inevitable that the rocket may fail to follow 
its prescribed path in space Errors in initial conditions or 
unknown gravitational influences will cause deviations from the 
predicted trajectory. 

As soon as differences between actual position(as described 
in Part A above and predicted position become known, 
a new situation arises. The decision must be made to either 
return to the previously calculated orbit by means of corrective 
thrust, or to recompute a new path which will be within the 
capabilities of the rocket. Time is an important consideration in 
this discussion since corrective action must place the rocket on 
a path which will result in almost simultaneous coincidence of 
Mars and rocket. This is discussed further below. 

It is presumed that almost continuous position finding data 
will be available so that comparisons can, in fact, be made 
between actual and predicted position. It is also reasonable to 
assume that deviations, hence corrections, will be incremental 
in nature. For the sake of simplicity, let us also assume that 
computation of a new trajectory, when required, can take place 
instantaneously. Under actual conditions, solution time can be 
allowed for with no great difficulty . 

The desirability of placing the rocket on a new path rather 
than forcing it to adhere to a previously defined path has been 
mentioned. The latter method implies, in effect, the application 
of continuous corrective thrust. Reference 26 has indicated an 



125 



analytical approach to the thrust requirements necessary to fly 
a prescribed path. Until such time as fuel supply poses no 
problem, however, continuous corrective thrust is prohibitive. 

As little application of corrective thrust as possible is desired, 
since there is little available for expenditure^ 

The objective of corrective action is not necessarily to 
achieve a pinpoint interception of Mars, but rather to reach a posi- 
tion favorable for entry into orbit about Mars. This might 
conceivably be a position at some distance thousands of miles 
from the planet as well as on the planet itself, since the 
gravitational influence of the arrival planet and decelerating 
thrust will combine to make entry into orbit a less significant 
consideration. Therefore, some error can be tolerated. However, 
when and if deviations occur which are outside tolerance limit, 
corrective thrust becomes necessary. Obviously, the tolerable 
deviation is smaller at the early stages of flight. 

For purposes of analysis let us assume that the allowable 
tolerance has been exceeded and that corrective action is called 
for. The rocket is moving primarily under the attraction of a 
single body, the sun, which is justifiable in accordance with the 
calculations of section 3.6. We will also consider Mars as a 
massless point to be "intercepted" (in the sense defined above) 
at a specified time. 

28 

The subsequent analysis parallels that of Lawden 
although different in many respects. 

7.11 Analytical Development 

I 

Let r^ , 6 ^ define present position in space at time t^ . Let 
rj^ , define the position the rocket would occupy at time t^ had 
flown the correct path, , ©2 define the position of Mars at time 
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tg, the planned time of interception It is desired to intercept 
Mars at time t 2 via a route which obviates the necessity of re- 
turning to the original trajectory. Figure 7,7 is helpful in the 
interpretation to follow. The route is a conic passing through 
P and J , 




Fig. 7, 7 Interception at Originally Predicted 

Time, t 2 

Clearly, it is a restriction to be forced to intercept the 
planet at time, t 2 , but results in a great simplification in the 
computer requirements. Let us pursue this analysis first. 

The conic path to be flown is presumed to be a segment of 
an- ellipse (because of velocity requirements) with sun as focus. 
The line of apsides of this new ellipse, whose elements are 
designated (a 2 j ^ 2 ) is aligned at an angle u with the reference 
line of the originally prescribed ellipse (elements apC^). 

From the fundamental elliptical equation 4-1 

a^ (1 - 

r = ? S (7-11) 

1 + e^ cos (0 - w ) 

^ o 
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(7-12) 
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2 



f 



^2 

1 + 02 



(02 ~ ^ 



f 



where r , r„, 0 , and 0„ are known. ^ 

Position in an elliptical orbit at any time t can be found from 
47 

the relationship 



/ 2 ' 
2a~t ev 1 - e sin 0j 

T i ‘ 

1 + Oj cos 0J 



I 

Cj + cos 0J^ 



t 

1 + e cos 0j 




where t is measured from t = 0 which occurs at perihelion passage 
of ellipse being considered.. 
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T, the period of the ellipse^can be expressed as a function 

of a, and is usually known. 

3/2 



2 



T = 
2 



\fcM 

From these relationships. 
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sin 



j + cos(0 2 ” w ) 
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l+Cg cos(02 - u ) 



It 

2 



1 + e„ cos (0 ~ u ) 

^ o 



, e„+ cos (0 “ u) 

. ~i 2 o , Tt 

sin +~2 

1+02 cos (02 - w) 



where tr, and t are unknowns, measured from the reference line 
2 o 

inclined at an unknown angle u to a known reference line. 



K the time of travel from r to r^ is constrained to be the 

o 2 

originally planned interval of time from r^^ to r 2 > the interval 
t 2 - t^ becomes a known quantity, and t 2 and t^ need not be known 
individually . 
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tr. - t = At - 
2 o 



62n/i - ^2^ sin (02 ^ ^o” 

i *“ i 

1 + 02 cos ( 02 “ w) ^ ^2 w) 

f I 

“1 e^ + cos (&2 ~io) e2+cos(0^- u) 

“ sin j 4sin 1 

1+ 62 cos (02 “ tj) 1+ e2cos(0^-u) 

, (7-14) 

In summary, position finding data will provide r^, 0^. The 
quantities r , 0^^ and r2^ 02 a-^'e known in advance as functions of 
time. Thus, At, the time till interception, is known for any r, , 

I " - 

9 ^. Equations 7 11 7-14 involve four unknowns a2, ^2, w, and T2, 

and are solvable by iterative methods . 

■* 

Since corrections are assumed to be incremental in nature 
when compared with the elements of the path being flown(a^, e^, 
^o^ ' ^2" ^2 replaced by a + 6a, e + 6e, and 

CO + 6(j. It is assumed T2 has been replaced by its value in 
terms of a2 . 

If these substitutions are made and small angle assumptions 
are valid (sin 610 = 6co, cos 6to= 1) equations in the three unknown 
6a, 6e, and 610 result, a^, e^ and co^ (assumed zero, initially) 
can be found as described in the navigational section and are good 
first approximations to a , e and co . 

Successive improvements of the accuracy of 6a, 6e, and 6to 
can be made if desired 

7.12 Corrective Action 

Since changes must inevitably be reflected in an alteration of 
velocity and angle of fire, we may again write two basic re- 
lationships 




(v^ + Av)^ 
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tan {6 + Ad) 



I 

sin ( 02 " ^ 2 ^ 

f 

1 - 02 cos ($2 - ^ 2 ^ 



Replacing a 2 , ©2 and 1^2 by a+ 6 a, 62 + 6 e and w+ 6 w, we 
arrive at 



(v^ + Av)' 



= GM [ f- 
^2 



_J 

a + 6 a 



Av 




1 

a + 6 a 




where is the velocity at the time corrective action is deemed 
necessary, and is the corresponding angle of orientation as 
defined earlier. Similarly, direction must be altered to effect 
A^^ which can be written 



(e +6e) sin (0o " w - 6w) 

A</> = tan ■ i 6 

l-(e+6e) cos (02“W“6w) 



a + 6a, e + 6e, and 10 + 6w are presumed to be known accurately 
as described in the preceding paragraphs. 

Corrective thrust must be actuated to produce the changes 
required. 



7.13 Interception at Newly Calculated Time 

Let us now consider the prospect of intercepting Mars at some 
time other than the originally planned time. 
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Mars 

Path 




Fig. 7.8 Interception at Newly Calculated 

Time 



It IS recognized at t^ that present position, r^, is un- 
satisfactory, i. e, a continuation along present track will not 
achieve interception. Moreover, assume that rocket corrective 
thrust capabilities will not permit an interception at r2> 0*2^ the 
originally planned arrival point. We desire a conic (elliptical) 
path from r^, 0 ’^ to some arbitrarily selected Mars position, r^, 
0*2 which has been chosen because it is within the rocket's thrust 
capabilities This new position may also, in fact, represent the 
position of Mars which can be most economically reached. The 
choice of this position does present a formidable problem, since 
the optimization process is necessarily involved. The selection 
of a Mars position is not discussed here although its importance 
IS emphasized. 



Presumably, there is on board a means of establishing Mars 

position as a function of some base time = a Mars Ephemeris of 

sorts. Fixing the desired interception point will fix the time, t^, 

of interception, t^ is known from an accurate time piece. Thus, 

the interval of corrective travel becomes do “ t . These times 

3 o 

are both based upon a common reference. The time interval is 
therefore determinable. 

Employing the previously developed reasoning procedures, 

j, = i 

1 + e„ cos {6 +u) 
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1 + e„ cos (6 + u ) 
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-1 ®2^ ^ -1 ^2^ ^^o ^ ^ ^ 

- sin j fsin i 



1+02 cos ( 0 g + w) 



1+ e„ cos (9 +u ) 
^ o 



and 



2w a, 



3/2 






By the methods previously described, the elements of the 
new path a 2 , 62 , W 2 can be found and corresponding 6 a, 6 e, and 
6 w's determined. It is now a relatively simple matter to 
calculate the velocity and angular changes necessary to fly the 
corrected path. 
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7 . 14 Summary 



It is seen that both approaches are essentially the same. 
From the standpoint of instrumentation, however, interception 
at time other than planned requires on board a knowledge of Mars 
position at all times and necessitates the selection of a new inter- 
ception point. The choice of this point is essentially an 
optimization process. 



CHAPTER 8 



GUIDANCE AND CONTROL 



8. 1 General 

The instrumentation of any guidance and control system for 
interplanetary travel is extremely complex. As stated previ- 
ously, this paper studies the so-called "combination" trajectory 
(i.e. - an initial thrust period followed by unpowered free-fall 
flight; modified, if neoessary, by corrective thrusts). Primary 
emphasis has been placed upon what has been defined as the "mid- 
course" of this type trajectory. Although considerable thought 
has been devoted to the use of departure and arrival orbits, the 
applicability of this paper is not limited solely to interorbital 
trajectories . Whether the vehicle orbits prior to commencing 
mid-course travel or simply continues directly from the ascent 
from the surface of the Earth does not affect requirements of the 
vehicle velocity vector orientation and magnitude and other para- 
meters for the actual transfer period. 

This Chapter outlines in simple form a suggested Mid- 
course Guidance and Control system based upon the material 
developed in previous portions of this paper. The major problem 
of reducing the system to the least possible weight and size has 
been ignored except to the extent that simplification has been 
attempted where it will not detrimentally affect the overall ac- 
complishment of the mission (i.e. - the arrival of the rocket at 
Mars within certain limits). 

Time-sharing of a miniaturized low -power digital computer, 
a common ephemeris memory, and one basic time generating 
mechanism have been assumed. This appears justifiable since 
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the elasped time between computations need not be of the order of 
fractions of a second, but can be several minutes or possibly hours. 
Digital computation is necessary due to the extremely low rates of 
change of rocket position, velocity, and acceleration; and the need 
for a high degree of accuracy. 

The problems of instrumentation have been divided into 
several sections as follows: 

(1) The Basic System 

(2) The Tracking Section 

(3) The Computing Section 

(a) Tracking Line Computation 

(b) Prediction Computation 

(c) Comparing Function 

(d) Correction Computation 

(4) Thrust and Stabilization Section 

(5) Errors and Effects 

These sections are treated below. 

8 . 2 The Basic System 

The functional responsibilities of Guidance and Control have 
been segmented into three basic categories, as shown in Fig. 8-1. 

A brief general description follows; 

The Tracking Section - Maintains Lines -of-Sight upon 

prescribed bodies and provides 
Tracking Line data therefrom to 
the Computing Section. 

The Computing Section - Performs the necessary computations 

to the Tracking Lines to obtain 
signals proportional to the cor- 
rections required to the rocket 
velocity vector to insure ac- 
complishment of the mission. 



COMPUTING SECTION 



•D 

C 

iP 

a* 

CO 




O) 

> 



o 

<u 

o. 

0 ) 



"O 

a> 



o 

Q> 



O 

o 



*o 

c 

o 



T 3 

<D 

> 



0 > 

(/> 

JD 

O 



T 3 

© 

<3 



*D 

© 



O 



© 

© 



o 

o 



</> 

X) 

o 



© 

k- 

o. 



© 

o. 

'u. 

o 

© 

X) 

3 

CO 



liJ 

I- 

o 



T 3 

a; 

w 

o 

a 

o 

a 

03 



S 

o3 

tuo 

o 3 

•fH 

O 

r-H 

C 

O 

•iH 

s 

o 

•rH 



03 

a 

o 



o 



(U 

SI! 



J 



a 

(V 

4-> 

w 

;>> 

w 



o 

Sh 



s:: 



o 

a 



T 3 

SH 

03 



(U 

a 

sH 

03 

T 3 



US 

tuo 



CO 

•iH 



137 



The Thrust and Stabi- 
lization Section - Provides vehicle orientation and 

corrective thrust in accordance 
with the received signals . 

The aforementioned system sections will be examined in 
more detail in the following sections of this paper. The Computing 
Section of the rocket involves such a complex problem that it has 
been subdivided into four parts, each of which will be discussed in 
its own Section of this Chapter. 

8. 3 The Tracking Section 

The tracking section envisioned determines tracking lines 
from the vehicle to the Sun and planets measured in an helio- 
centric ecliptic coordinate system. These angular measurements 
are then sent as electrical signals (TL^ to the tracking line com- 
puter in the tracking section. Two systems are presented which 
differ in method of space reference frame determination, number 
of star and planet trackers required, and in probable degrees of 
accuracy and reliability, 

8. 31 System One 

The first system obtains space reference from star trackers 
rigidly mounted to a control member platform isolated from vehi- 
cular motion by gimbals. This platform is "fixed" in space and is 
maintained in this attitude by tracking error signals generated in 
the fixed star tracking units, which are used to actuate the gimbal 
drive motors . A gyro package is mounted on the controlled member 
solely to assist in holding track during periods of radical maneuver 
and to aid in initial lock on. Three additional star trackers are 
mounted on the controlled member, each with its own tracking 
member drive system, as illustrated in Fig. 8,2 and 8.3. One 
tracking head determines the rocket-sun tracking line while the 
others determine two different rocket-planet tracking lines. 
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Fig. 8.2 



Illustrative schematic of tracking system 
with gyro assist 
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Fig. 8.3 Line functional information flow diagram of Tracking System I 




Angular measurements of these lines are made in heliocentric 
ecliptic coordinates, (see Fig, 7.1 ) and transmitted to two 

separate range computers within the tracking line computer of 
the computing section. 

8. 32 System Two 

The second tracking system (Figs. 8,4 and 8,5) consists of 
a space reference gyro package, and only three star trackers. 

The reduction in the number of trackers used is made solely to 
simplify construction at the expense of system accuracy Three 
orthogonally positioned single degree of freedom gyros are as- 
sisted by a single "fixed star" tracker rigidly mounted to the 
controlled member. Smoothed error signals from this tracker 
are combined with the signals from the gyros. This resultant 
signal IS then used to actuate the gimbal drive motors maintaining 
the controlled member reference plane in proper orientation. 

Two star trackers are then mounted on the controlled member, 
each with its own tracking member drive system. As in the first 
system one tracker determines the rocket -sun line, while the 
second tracker determines a rocket-planet line. Angular measure- 
ments of these lines feed into a range computation section of the 
tracking line computer, 

8. 33 Comparison of Systems One a nd Two 

It IS obvious from comparing Figs . 8.2, 8,3, 8.4 and 8.5, 
that the second method discussed in 8.32 is much simpler to 
design and construct since only three stellar tracking heads are 
needed while the system discussed in 8.31 requires six tracking 
heads (Possibly this could be reduced to five if only two heads 
v/ere used for space reference determination). The main draw 
back to system two is the dependence upon gyros for space refer- 
ence. These units are not, at present, stable enough nor pre- 
dictable to the degree required, for a time of flight of the length 
required for an Earth -Mars transit. Section 8,64 discusses this 



\l/ 




Fig. 8.4 Illustrative schematic diagram of tracking system 
with stellar assisted gyro space reference 
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Fig. 8. 5 Line functional information flow diagram. Tracking Section II 



particular aspect, 

A second drawback to the simplified system of section 8.32 
is the use of only one planet The effect of this is felt in com- 
putation of ^^^OA^obs" treated in^section 8,41. 

8 4 The Computing Section 
8 . 41 Tracking Line Computation 

The tracking line computation section of the Guidance System, 
receives the angular measurements of the tracking lines de- 
termined in the tracking section. Using these measurements and 
data from the Ephimeris Storage and Time Generator within the 
computing section;, ^I^QA^obs ^^©A^obs determined. In 
addition the various separate elements of the ellipse are com- 
puted for complete identification of the actual trajectory being 
followed by the vehicle. 

The method of computation is geometric using discrete 

measurements, and not the integration of velocities and acceler- 

71 72 

ations experienced by the vehicle. ' The reasons for this 
are shown in section 8.6. Chapter 7 developed the equations 
necessary for the solution of this problem and illustrated the 
various geometric relations in Figs. 7, 2 through 7. 6, 

The first tracking system discussed in section 8. 31 has as 
outputs TL0 , TL0 and TL(~^ , These are inputs to two separate 
range computers which determine two values of F-0^ using the 
equation of section7.6 . These two values are then averaged and 
smoothed for greater accuracy and then combined with the vehicular 
heliocentric latitude and longitude (L^, as presented by TL0 , 
Since a series of must be used in computation of the elements 

of the elliptical trajectory, these must be held in a memory storage 
unit, and fed to the trajectory elements computer upon demand. 

The determination of CV0^)^|^g is made by a simple integration 
process and along with ^ direct output of the tracking 

line computer. The output of the trajectory elements computer 
are 0, a, e, and i observed. The simplified iriformation flow 
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diagram of the above system is shown in Fig. 8,6 

The second system is the same as that described above but 
uses only one range computer with the resultant loss of a degree 
of accuracy in the solution, although a system simplification is 
achieved which may be permissible. A simplified information 
flow diagram of this is illustrated in Fig. 8, 7 



8 . 42 The Prediction Computation Section 

The Prediction Computation Section consists primarily of a 
number of components collectively termed the Prediction Com- 
puter, (see Fig. 8,8). The purpose of the Prediction Computer 
is to continuously establish, as a function of time after t^, the 
position and velocity of the rocket relative to the reference body. 
The position and velocity so established may be termed "pre- 
dicted" in the sense that they do not stem directly from obser- 
vations, but rather are values which should actually occur if no 
uncompensated events transpire to alter the rocket motion. The 
computer must establish these values without reference to any 
source of information from outside the rocket itself. This in no 
way prohibits the computer from sharing storage data, power 
supply, time generation, and the like, with other components. 

If no corrections to the trajectory were envisaged all the outputs 
could be computed prior to launch, stored in a memory, and 
sampled as needed, Since corrections are envisaged the computer 
must possess the capability of solving values for the rocket po- 
sition and velocity during the entire time of flight. 

Selecting the Sun as the center of our inertial reference 
frame we may write Equation 3-11, 



^OA 



j = ^nf ^©A ■ %q) + a " ^<Iq) + ^^dA ' 



(8-1) 
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':^ig. 8.6 Line functional information flow diagram for tracking line computer 

(System One) 
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Fig. 8. 7 Line functional information flow diagram for tracking line computer 

(System Two) 
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From Fig 2.8 we determine 



G^O < 10"^ ft/sec^ 



G^0 < 10 ft/sec^ 





= 0 



(8-2) 



Hence the equation of motion to be solved continuously by the 
Prediction Computer is 



The threshold sensitivities of the components will determine 
the Phases referred to in Section 3 6. Below these thresholds 
a particular component will be inactive and the problem will be 
reduced to the appropriate one of Equations 3-14, 3-15, or 3-16. 

Referring to Fig 8 8 the functions of the sub -components 
are almost self-explanatory when examined with regard to 
Equation 8 '3. The Summer forms the quantity representing the 
acceleration of the rocket with respect to the Sun, The acceler- 
ation is integrated twice to form the rocket position vector, 
which is fed-ba,ck to obtain values for those gravity fields con- 
sidered in the equation. Thrust and rocket mass form the 
remaining inputs tothe Summ.er 

43 Comparator Section 

The comparator section may be considered an integral part 
of the main computer Its function is to compare observed 
position in space, emanating from the tracking 

computer, with predicted position in space, derived 

from the prediction computer. Deviations of observed position 



«OAjy fnf+G0A + °®A"GaA^GcfA 



(8-3) 



from predicted position are de' 




suiiicient 



magnitude, form the basis for corrective action. The test of 
"sufficient magnitude" is a function of the position of the vehicle 
with respect to Mars. Obviously, the permissible deviations 
would be small during the early stages of flight, when the position 
of the rocket relative to Mars is quite lar'ge . Similarly, as this 
relative position becomes less, tolerances become larger and 
errors assume a lesser importance. A deviation computer 
weighs or assesses the rocket position relative to Mars and 
creates this standard of measurement. Fig. 8.9 is a functional 
representation of this section. 

If actual deviations in position, (D^) exceed those 

permissible deviations, (D ) RQa, the correction comparator 
is assumed to emit a signal actuating the corrector section. 

This will automatically cause initiation of the computation of 
the new trajectory. 

8. 44 The Corrector Section 

The corrector section can properly be considered a most 
complex and critical part of the computer. Upon its proper 
performance will depend an accurate approach to the objective 
planet or the alternative, a body in free orbit in space forever. 

The corrector section is continually supplied with observed 
rocket position, elements of the observed 

trajectory, as well as the position of Mars with respect to the 
Sun as afunction of time, RQ^t). These quantities form the 
basis for the computation of a new path to the sought-after planet. 
A discussion of this general problem has appeared in Chapter 7. 

The initiation of operations within this section will be de- 
pendent upon receipt of the comparator section actuating signal, 
Sg(Act). 

The modified trajectory is assumed chosen through an 
optimization process which minimizes fuel expenditure, time. 
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Fig. 8.9 Line functional information flow diagram for comparator section 
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Fig. 8,10 Line functional information flow diagram for corrector section 



or some predesigned parameter. Once this selection process 
is completed, signals proportional to the required velocity 
changes can be converted into control system commands, (S )T, 
which will eventually translate the changes into corrective 
action. (See Fig 8.10. ) 

From a practical viewpoint, however, the selection of an 
alternate path in flight can be expected to pose a difficult 
problem Size and weight considerations will necessitate an 
absolute economy in functional employment of components. Even 
then, the optimization problem is complex and may well repre- 
sent a major obstacle. 

8 5 Thrust and Stabilization 

The thrust and stabilization unit is essentially a modification 
of the familiar control loop, altered because of operating en- 
vironment to include rocket thrust device rather than control 
surfaces, A signal proportional to the required velocity actuates 
the thrust device which is positioned prior to firing. Several 
means of directional control are available for corrective thrust 
requirements ; 

(1) Movable vanes in the jet stream 

(2) Small, individual variable thrust jets 

(3) Universally mounted main jet engine 

(4) Gimbal mounted rocket thrust chamber 

Posiuoning of the thrust device would be in response to the 
correct or section output signal, as would the controlled 
metering of propellant to insure a velocity change of proper 
magnitude. 

The accelerated motion of the rocket produced by the thrust 
device would then be sensed by an accelerometer which com~ 
pares, in effect, commanded thrust, with resulting thrust. 

The stabilization of the rocket in space is another problem 



of no little difficulty. Some of the proposals for rocket orien- 
tation have included micro-jet control, gyroscopic torque 
stabilization, and flywheels . Many authors combine vehicle 
motion and the stabilization although they distinctly differ and 
require functional separation. 

Fig. 8.11 is a line functional information flow diagram of a 
proposed thrust and stabilization sections. 

8.6 Errors and Effects 

It has been stated previously that the unpowered elliptical 
orbit trajectory for transfer between planets is subject to large 
errors in terminal position caused by very small perturbing 
forces during transit and also very slight errors in velocity 
vector orientation at the end of the initial thrust period. If, as 
has been assumed in this paper, corrective action must be taken 
during mid-course, the reduction in terminal error will be a 
function of time remaining till intercept, accuracy in determining 
the actual trajectory being traveled and the accuracy of determi- 
nation and orientation of the corrective velocity vector. We will 
now show the magnitude of some of the terminal errors caused 
by various assumptions and inaccuracies . 

8. 61 The Assumption of Circular Coplanar Orbits 

The assumption of coplanar circular orbits has been made 
by numerous authors for the sake of simplicity in presenting 
basic ideas. > a,nd others. If we compare the 

orbits of Earth and Mars and transfer between these orbits we 
find the following to be true for the worst possible cases. 

At "favorable opposition" (once every 16 years) Mars is 
6 

35 X 10 miles from Earth. At "unfavorable opposition" Mars is 
0 

63 X 10 miles from Earth. This gives a maximum actual vari- 

0 

ation in distance between orbits of 28 x lO" miles (. 3 AU), while 

0 

the so called average distance is only 49 x 10 miles (. 523 AU). 
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This is a 57% variation which is ignored with the assumption of 

circular orbits. Using the coplanar assumption (not taking the 

actual orbital inclination of 1° 51" into account) we introduce a 

0 

terminal vectical error of 4.57 x 10 miles for the worst possible 

V 

case. 

8. 62 Errors in Astronomical Distance 

The basic unit of astronomical distance is the Astonomical 

Unit (AU). This is determined by the mean range of the sun from 

the Earth and by international agreement (Paris, 1911) has been 

6 

established to be 93 x 10 miles . However, this distance is 

based on angular measurements of the diameter of the sun. Due 

to limits of resolution, limited base lengths for triangulation 

and other effects, this basic unit is thought to have a probable 

54 69 

error of 10,000 miles. This uncertainty is therefore 

carried over to the measurement of positions of solar bodies 
and hence to vehicular orientation in solar space. Until more 
accurate astronomical measurements are available a rather large 
residual error in spatial position will have to be accepted and 
worked with. 

8. 63 Errors in Initial Velocity Vector 

If we now take a typical trajectory from Earth to Mars, 
we can see the effects of errors in initial velocity vector on 

0 

terminal positions. Assume a trajectory of arc length 354 x 10° 
miles such as used in section 4.14 . This trajectory has a time 
of flight of 172 days, and an orbital inclination of 2. 51°. An initial 
velocity error of 1 foot per second for this trajectory will give 
an error in arc length alone of 2820 miles. This is neglecting 
the effects on the actual elliptical properties of the orbit which 
increase the actual miss distance even more. Using an average 
orbital velocity of 23. 8 miles per second, this would be equiva- 
lent to an error in arrival of 118.5 seconds, at which time Mars 
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(traveling approximately 15 miles per second) would be 1780 miles 
behind the initially proposed intercept point. 

An error in inclination or angular orientation of the initial 
velocity vector of . 01 degrees gives a terminal error normal to 
the trajectory of 18^ 600 miles . We can also approach initial 
angular error from the standpoint of what will be the requirement 
for insuring a terminal error of not more than 1, 000 miles. This 
will necessitate an accuracy within 1. 94 second of arc. 

49 

J. G. Porter discusses the effects of small initial velocity 
vector errors on the shape of a simple Earth-Mars elliptical 
trajectory. This trajectory starts at its perihelion and arrives 
at Mars at aphelion and has a basic time of flight of 281 days. Porter 
concludes that in order for this trajectory to arrive at Mars 
within a distance of 50, 000 miles, the velocity must be controlled 
to .001 miles per second or one part in 18, 000. This would mean 
that to get within 1, 000 miles of the desired terminal position the 
initial error would have to be less than .106 feet per second. 

For this same trajectory Porter concludes that an error 
of 12 seconds of arc in initial velocity vector orientation will lead 
to a 50, 000 mile error at termination. This means that for a 
1, 000 mile error we would have to be within . 4 of a second of 
arc at burn out. 

The major conclusion which can be drawn from a look at 
errors in initial velocity vector orientation is that a simple 
unpowered trajectory with no corrective thrust applied during 
mid-course will not suffice for interplanetary travel. This was 
previously stated in section, 4. 4 HDwever, of the two errors con- 

sidered, velocity and angular, the errors in angular orientation 
are the more critical. 

8 , 64 Mid-Course Errors 



We can now look at errors in some instruments which might 



be used within the vehicle for guidance purposes such as star 
trackers, accelerometers, gyro packages, and associated 
equipment. 

If we presume the need for a maximum star tracker error 
equal to or less than 1, 000 miles in determining the positions of 
solar bodies 1 AU distant, we obtain a limit of 2.22 seconds or 
.0382 minutes of arc. This is a maximum allowable closed loop 
error composed of errors in reference plane, tracking head, 
and all associated equipment. 

Another approach to this problem of star tracker error is 
to assume that each sighting is a random event with Gaussian 
distribution of the error in angular measurement. K we now 
postulate that we have a "maximum" error of -5 seconds of arc 
in the values obtained from a star tracker and that 99.7% of the 
measurements will lie within these values, 95% will be in error 
less than 3,33 seconds while 68% will be in error less than 1.67 
seconds. Since an error of 1 second of arc is equivalent to 450 
miles at a distance of 1 AU, we can assume that 68% of the 
measurements would have an error less than 750 miles, and 95% 
would have an error less than 1, 500 miles when sighting bodies 
1 AU distant. It appears, then, if a number of readings are 
taken over a short period of time with proper smoothing and 
averaging, suitably accurate values can be obtained from tracker 
systems which may by themselves have a considerable possible 
error. 

To determine an upper limit of the threshold of sensitivity 

to acceleration required for the trajectory of Section 4.14 of this 

paper we can again use a maximum error of 1, 000 miles at 

-»8 

termination. This corresponds to 4.75 x 10 feet per second 

for a time of flight of 172 days. For this same trajectory the 

maximum allowable drift of a purely "inertial" guidance system 

to insure vertical displacement error of not more than 1, 000 miles 

~5 

would be 1. 64 X 10 mr per day. This is determined using the 
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0 

trajectory arc length of 354 x 10 miles and time of flight of 172 
days . 

The primary conclusion that can be drawn from the above 
simplified mid-course guidance system errors is that although 
completely self contained inertial guidance systems are highly 
desirable^ the required accuracies are rather high; particularly 
if return trips are contemplated without rezeroing the system. 
However, with the addition of star tracking units, suitably 
accurate guidance can be obtained which is not limited by or 
dependent upon time of flight. 
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CHAPTER 9 



CONCLUSIONS AND RECOMMENDATIONS 
FOR FURTHER STUDY 



9.1 Introduction 

The contents of this paper lead to several conclusions and 
recommendations of interest to those pursuing the problem of 
interplanetary travel. Many areas included in this thesis are 
presented with a view toward future use by students of the field 
and lead to a better definition of the problem. The bibliography 
represents a comprehensive literature search of available un- 
classified related material. 

9, 2 Conclusions 

(1) During the time of flight of a rocket from Earth to 
Mars the center of the Sun may be considered an 
inertial point. 

(2) During a major portion of the trajectory the problem 
may be considered one of two bodies^ with some 
perturbations . 

(3) Threshold sensitivities and drift uncertainties of 
current gyroscopic mechanisms preclude the use of 
a pure inertial navigation system at this time, par- 
ticularly when a return trip is contemplated. 

(4) During the time of flight of the rocket the motion of the 
"fixed" stars is so small that they provide a means of 
establishing a space reference frame. 

(5) Determination of rocket present position in space within 
the limits proposed in this paper lirs within the capa- 



161 



bilities of the current "state of the art". 

(6) Within the limits of the assumptions made, the 
"combination" trajectory fi p ,an initial thrust period 
followed by free fall flight; modified, if necessary, by 
corrective thrust) appears to offer the most immedi- 
ately feasible method of travel to Mars . 

(7) Rocket velocity vector orientation appears to be more 
critical than magnitude. 

22 

(8) The Hohmann ellipse provides the minimum required 

velocities to get to Mars but if rocket capabilities 
include (as defined within this paper) of 10, 000 ft /sec 

or more a wide range of trajectories is available. 

(9) The Rocket Plane Concept provides a means of reducing 
the trajectory to a two-dimensional problem, thereby 
greatly simplifying rocket construction requirements. 

It is emphasized that certain assumptions often made 
in the literature to reduce the problem to two -dimensions 
are unrealistic when viewed in the light of actual 
condition. The Rocket plane concept does not assume 
two -dimensions but selects the reference plane such 
that deviations from this reference plane are negligibly 
small. 

(10) On a kinetic energy basis a certain optimum departure 
orbit exists for any given trajectory to Mars. Similarly 
an optimum arrival orbit exists about Mars , Additionally 
these optimum orbit radii are characterized by the 
property that for a relatively small increase in available 
velocity a satisfactory departure orbit radius may be 
used which is much nearer to the Earth surface than 

the optimum orbit radius, 

(11) Instrumentation for travel between planets appears to 

lie within the capabilities of the present "state of the art". 
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9. 3 Recommendations for Further Study 

(1) Digital computer trajectory analysis should be 
instituted covering: 

(a) Required parameters for all trajectories whose 
initial rocket velocity lies within reasonable limits. 

(b) The effect of variations in initial conditions on 
linear miss distance at Mars. 

(c) Perturbative effects of solar system masses which 
were ignored in this paper due to their small 
magnitude. 

(2) Additional study should be instituted on the design and 
construction of free-fall space reference stabilization 
systems suitable for interplanetary travel. 

(3) Studies should be commenced which would lead to the 
establishment of optimization standards for trajectories 
to Mars on the basis of both fuel expenditure and 
kinetic energy. 

(4) A study of permissible deviations from prescribed path 
during flight allowing terminal accuracy within pre" 
scribed limits should be instituted. 
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SYMBOLS 

DEFINITION SUMMARY 1 



Symbol 




Definition 


© 




Sun 


c 




Moon 


0 




Earth 


cf 




Mars 


z 




Jupiter 


\ 




Saturn 


6 




Uranus 






Nep tune 


E 




Pluto 


? 




Mercury 


9 




Venus 


☆ 




Star 



a 


Semi major axis of any ellipse 


b 


Semi minor axis of any ellipse 


c 


The magnitude of the effective 
exhaust velocity vector measured 
with respect to the walls of the 
rocket nozzle 


e 


Eccentricity of any ellipse: 



defined by the relation 




a 

Unit field forces - forces per 
unit mass whose existence is due 
to the presence of a field 
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Symbol 



f 



nf 



i 



j 



m 



A 



P 



r . 
D 




t 

t 



b 



t 

c 



Definition 

Unit non field force - forces per 
unit mass whose existence is due 
to all causes other than fields 

Inclination of orbit - the angle 
measured from reference plane 
to plane of orbit, in a plane 
perpendicular to the line of nodes; 
generally positive where measured 
in a direction "above" the refer- 
ence plane 

A variable indicator which may 
take on identities as propounded 
within the text 

Instantaneous mass of rocket 

An arbitrary planet 

Magnitude of radius vector from 
the reference focus of an ellipse 
to the indicated point on that 
ellipse. Subscript j refers to 
indicated quantity 

Radius of the circular orbit about 
the indicated body 

The independent variable time 

Burning time - that period of 
time wherein a given amount of 
propellant is consumed 

The time at which an imaginary 
rocket would have arrived at Mars 
future position if it had departed 
aphelion of its ellipse at t = t^ 

Note; t- = t - t , 
fed 
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Symbol 



Definition 



t 

o 

A 

AU 



D 

E 



e 

F 



The time at which an imaginary 
rocket would have arrived at 
Earth present position if it had 
departed aphelion of its ellipse 
at t = t 

o 

Time of flight 

Time at "j" instant; subscript 
"j" refers to indicated quantity 

The instant of time at which a 
given operation commences . 

Angle 

Astronomical Unit - the distance 
associated with the mean parallax 
of the sun. By international a- 
greement (Paris 1911), this was 
set as 8.80" of arc and corre- 
sponds to about 93, 000, 000 miles. 
In this text the numerical value of 
the AU will be taken as 93, 000, 000 
miles . 

Drag - as conventionally defined 

The ecliptic plane - the plane 
defined by the loci of points 
representing the instantaneous 
position of the earth throughout 
the year measured with respect 
to the sun 

Eccentric Anomaly 

Force - as conventionally defined 



lo7 



Symbol 

G 






GM. 

] 



Imp 



KE 

L. 

] 



L 

M 

P 

M 

o 



Definition 

Gaussian gravitational constant 

approximately equal to; 

-8 3 2 

6.67 X 10 " cm /gram (m) sec 

Directed acceleration due to 
gravity field of the indicated 
body. A second subscript refers 
to that body experiencing this 
directed acceleration. Subscript 
"j" refers to indicated body 

The product of the Gaussian 
gravitational constant and the 
mass of the indicated body. Sub- 
script "j" refers to indicated 
body 

Impulse - the integral of the 
magnitude of thrust over the 
time of its duration 

Kinetic Energy - as conventionally 
defined 

Celestial heliocentric latitude - 
angle measured in a plane perpen- 
dicular to the ecliptic plane from 
the ecliptic plane to the line drawn 
from the sun to the indicated point; 
measured positively toward ce- 
lestial north 

Lift - as conventionally defined 

Total mass of propellant 

Gross mass of rocket at instant 
of firing 
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Symbol 

NN' 

PE 

PR 

R 

R 

"t 

TE 







Definition 

Line of nodes - line drawn through 
the ascending and descending nodes 

Potential Energy - as conventionally 
defined 

Payload Ratio - PR = 1 - ^ 

Vector from the reference point 
to the indicated point 

Magnitude of R 

Thrust - reaction force on the 
rocket or vehicle unit structure 

Total Energy - as conventionally 
defined 

Velocity vector of the indicated 
body with respect to a heliocentric 
non rotating coordinate system. 

Velocity vector of the rocket 
relative to the indicated body, 
and measured in a heliocentric 
non rotating coordinate system 

Velocity vector of the rocket in 
circular orbit about the indicated 
body measured in a body centered 
non rotating coordinate system 

Velocity vector of "escape" from 
the gravity field of the indicated 
body, equal to 

'sfT X T. 

3S 

measured in a body centered non 
rotating coordinate system 
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Symbol 



Definition 



■^(j A)g 



V. 

me 



Velocity vector of the rocket with 
respect to the indicated body 
taking into account the kinetic 
energy required to'' overcome the 
effect of the gravity field of the 
indicated body^ measured in 
heliocentric non rotating coordi- 
nate system 



That increment of V/. which 
must be provided by the rocket 



a 

7 



I 



e 



Angle measured from the x-axis 
to the thrust vector 

Angle measured from the y-axis 
to the thrust vector 

Angle measured from the z-axis 
to the thrust vector 

X, y, and z axes represent an 

orthogonal coordinate set. The 

thrust vector has its origin at 

the origin of this coordinate set 

M 

Propellant loading 'factor = 

Mq 

0-w for the orbit of indicated 
body 

The angle, measured in the rocket 
plane of motion, between a body- 
sun line and rocket ellipse aphelion- 
sun line; positive angles when di- 
rection of measurement is opposite 
to rotation of planets about the sun. 



170 



Symbol 



Definition 



Q (contd) 



6. 

] 



c 

e 



s 

X. 

3 



7T. 

3 



%A 



Prime (*) indicates measurement 
from perihelion vice aphelion with 
direction of measurement positive 
when same as rotation of planets 
about sun 

Q at any time j 

0 at t = t j 
d 

0 at t = t 

c 

0 of Mars at t = t , 
d 

Celestial heliocentric longitude 
- angle measured in the ecliptic 
plane from the reference direction 
Aries ( HT ), in the direction of 
rotation of planets about the sun, 
to the projection on the ecliptic 
plane of the line drawn from the 
sun to the indicated point. 

Longitude of perihelion of in- 
dicated orbit 

The angle^ measured in the 
rocket plane of motion included 
between the velocity vector of 
the rocket or vehicle and a line 
perpendicular to the radius vector 
from the sun to the vehicle 

<j> at t = 

<i at t = t 
^ c 

The angle measured in the rocket 
plane of motion included between 
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Symbol 



Definition 



%A 

■‘tJ'A 



OJ . 

3 



A 



rri 



T 



n. 

3 



V0A and the projection of Vq 0 
on the rocket plane. 

The angle measured in the rocket 
plane of motion, included between 
and the projection of 
on the rocket plane . 

The angle included between 
perihelion and the line of nodes 
measured in direction of ro- 
tation of planets for indicated 
orbit. 

Rocket or vehicle 

The period - elapsed time between 
repetitive events 

The vernal Equinox (the first point 
of Aries) defined as the direction 
indicated when projection of the 
line of intersection of the earth 
equatorial plane and the earth 
ecliptic plane passes through 
the center of the sun; this action 
occuring as the sun apparently 
passes from south to north of the 
equator 

Longitude of ascending node of 
indicated plane 

Longitude of descending node of 
indicated plane 

Ascending node - celestial heli- 
ocentric longitude occupied by a 
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Symbol 



Definition 



body passing from "below" to 
"above" the ecliptic plane. 
Descending node - Ascending 
node plus 180° "Above" ecliptic 
defined as that part of space 
wherein an observer would note 
counter-clockwise rotation of 
planet. 
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Earth angular position measured from Earth aphelion 
as a function of time (Overlay #1) 
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